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PREFACE 


Program  DSSMATH  enables  mathematicai  manipulation  of  data  stored  in  the 
Hydrologic  Engineering  Center’s  Data  Storage  System  (HEC-DSS).  The  program  provides 
capabilities  for:  arithmetic  computations,  transformations,  such  as  stage  to  flow;  screening; 
and  estimation  of  missing  or  erroneous  values.  The  program  may  be  used  in  an  automated 
batch  environment  for  processing  a  real-time  data  stream,  or  It  can  be  used  interactively  to 
perform  ad  hoc  operations. 

DSSMATH,  previously  known  as  WCCOMP,  was  developed  in  1986  for  the 
Sacramento  District's  Reservoir  Control  Section.  WCCOMP  was  developed  as  a  team  effort 
between  HEC  personnel  (Dennis  Huff,  Alfredo  Montalvo  and  Jeff  Houghten)  and  Sacramento 
personnel  (Paul  Pugner  and  Richard  Jerome).  Its  primary  function  in  1986  was  to  provide  the 
Sacramento  District  with  the  capability  of  doing  mathematical  manipulations  on  HEC-DSS  time 
series  data  related  to  the  operation  of  reservoir  control  systems.  In  October,  1991  WCCOMP 
was  renamed  DSSMATH  and  its  scope  of  application  broadened  to  encompass  a  wider  range 
of  mathematical  manipulations  of  HEC-DSS  data. 
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1.  Purpose 

DSSMATH  enables  mathematicai  manipulation  of  data  stored  in  the  Hydrologic 
Engineering  Center’s  Data  Storage  System  (HEC-OSS).  The  program  provides  capabilities 
for:  arithmetic  computations,  transfbrmaUorra,  such  as  stage  to  flow;  screening;  and 
estimation  of  missing  or  erroneous  values.  The  program  may  be  used  in  an  automated  batdi 
environment  for  processing  a  reaMime  data  stream,  or  it  can  be  used  interactively  to  perform 
ad  hoc  qaerations. 


2.  Description 

The  present  version  of  DSSMATH  is  essentiaily  identical  to  a  predecessor  program 
WCCOMP.  DSSMATH  also  encompasses  selected  routines  from  the  program  MATHPK.  This 
user’s  manual  is  valid  fbr  all  the  different  implementations  (Harris,  UNIX,  and  DOS). 

DSSMATH  supports  regular,  irregular-lnforval  HEC-DSS  time-series  data  and  HEC-DSS 
paired  function  data.  The  standard  memory  cortfiguration  (Harris,  UNIX,  and  DOS  Lahey 
Extended)  can  accommodate  10  time  series  records  of  up  to  5300  values  each,  5  paired-data 
records  each  with  30  dependent  variables  with  a  maximum  of  5000  values  per  record,  and  15 
scalars.  Data  may  be  accessed  in  several  HEC-DSS  files  at  the  same  time.  In  addition  to  the 
DOS  Lahey  Extended  implementation,  a  cut  down  DOS  version  that  runs  under  640K  is 
available,  its  configuration  is  5  time-series  records  consisting  of  3180  values  each  and  3 
paired  functions  consisting  of  5  sets  of  dependent  variables,  with  250  values  per  paired 
function. 

Time-series  records  and  look-up  t^ies  are  referenced  in  the  various  processing 
functions  by  user-assigned  labels.  The  labels  are  one  to  six  characters  long  and  are 
assigned  when  the  data  is  computed  (CO  command)  or  retrieved  (GET  command)  from  a 
HEC-DSS  file.  The  HEC-DSS  file  name  and  data  pathname  are  specified  as  parameters  to 
the  GET  commarKi  (See  Figure  1). 

In  either  an  interactive  or  in  a  batch  setting,  DSSMATH  proceeds  according  to  a 
sequence  of  instructions  provided  by  the  user.  A  typical  sequence  of  instructions  consists  of 
a  retrieval  of  data  from  the  HEC-DSS  file,  a  series  of  computations  involving  the  data,  and 
storage  of  results  in  a  HEC-DSS  file. 

Computations  in  DSSMATH  consist  of  arithmetic  manipulations  of  tinrte  series  data, 
such  as  multiplying  a  time  series  by  a  constant  or  another  time  series,  and  special  functions 
that  pro^^  more  complex  and  specific  capabilities,  such  as  rating  table  transformations. 

Elementary  screening  functions  flag  data  values  internally  to  fadtitate  subsequent 
corrections.  Replacements  fbr  flagged  and  missing  values  are  provided  by  simple,  linear 
estimation  functions.  Internal  flags  are  used  to  assure  that  missing  data  do  not  erroneously 
affect  computations.  However,  internal  flags  can  not  be  read  or  written  to  HEC-DSS  at 
present,  and  are  not  compatible  with  HEC-DSS  Version  6  flag  scheme. 
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The  plotting  capabilities  of  computer  program  DSPLAY  have  been  incorporated  into  a 
spedai  verdon  of  DSSMATH  named  D^MATHO.  These  include  control  of  the  plotting 
parameters  and  options,  such  as  grids  and  shading.  The  user  assigned  labeis  can  be  plotted 
by  DSPLAY  in  addition  to  DSPLAY’s  standard  displays  of  DSS  data.  The  majority  of  the 
DSPLAY  commands  are  available  with  the  exception  of  commands  that  would  conflict  with 
DSSMATH  commands.  DSSMATHD  is  only  available  on  the  Harris  Computer  and  at  the 
present  time  there  are  no  plans  that  this  capability  will  be  added  to  the  UNIX  aixi  DOS 
versions.  (Ref:  HEC-4DSS.  Users  Guide  a^  Utilitv  Program  Manuals^ 

The  features  of  PREAD  are  provided  in  DSSMATH.  PREAD  facilitates  efficient 
repetition  of  routine  procedures.  For  example,  a  complex  series  of  computations  can  be 
stored  in  a  PREAD  macro  arKi  started,  *or  example,  as  a  selection  from  a  screen  menu  or  a 
graphics  tablet  menu.  (Ref:  PREAD.  Functions.  Macros.  Menus  and  Screens  User 
Information^ 


3.  Use 

3.1  Running  the  Program 


DSSMATH  is  initiated  with:  DSSMATH  [parameter] ... 

"parameter*  has  the  form  "key  word^xxx"  and  may  be  one  or  more  of  the  following 


specifications: 

Keyword 

Default 

Descriotion 

INPUT 

*0 

Terminal  or  batch  job  input  file 

OUTPUT 

*3 

Terminal  or  batch  job  output  file 

FUNFILE 

MATHFUN 

PREAD  function  file 

MACFILE 

MATHMAC 

PREAD  macro  file 

LOGFILE 

W3 

PREAD  log  file 

SCNFILE 

GENSCN 

PREAD  screen  definition  file 

In  the  interactive  environment,  the  default  output  is  the  user's  terminal  display.  In  the 
batch  environment,  the  default  output  is  the  job  output  file. 


3.2  Use  of  Commands 

The  processing  of  timaeeries  data  with  DSSMATH,  either  interactively  or  through  a 
batch  input  file,  is  directed  with  commands.  Commands  are  specified  with  two-character 
mnemonics  in  columns’  1-2  of  the  input  stream.  The  generalized  form  of  the  commands  is: 
XX[.options]  [parameter  1]  [parameter  2] .... 

Commands  typically  require  parameters.  Parameters  supply  specific  information  about 
the  objects  or  entities  to  be  processed  or  specifications  to  control  the  process. 

Options  are  single  characters  used  to  specify  alternative  ways  of  processing. 
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Input  lines  beginning  with  two  asterisks  (**)  are  used  tor  notes  or  comments  and  are 
simply  echoed  in  the  output 

A  list  of  commands  is  presented  in  Table  1.  The  commands  are  described  in  detail  in 
the  appendices  A  and  B. 


Table  1  Commanda 

Command 

Description 

** 

Comment 

CATALOG 

Display  a  catalog  of  a  DSS  file. 

CLEAR 

Remove  a  variable  data  label  from  memory. 

COMPUTE 

Perform  a  computation. 

OIAG 

Toggles  DSS  diagnostic  trace  output  on  and  off. 

DPATH 

Display  a  selective  catalog  of  a  DSS  file. 

FINISH 

Terminate  and  exit  the  program. 

GET 

Retrieve  data  from  a  DSS  file. 

HELP 

List  commands. 

PLOT 

Plot  veri-^les. 

PUT 

Store  data  in  a  DSS  file. 

SO 

Set  data  descriptions. 

SHOW 

Displays  selected  internal  information  about  the  data  variable  denoted 
by  specified  variable  label. 

SMiSSING 

Set  missing  value  indicators. 

SP 

Set  data  pathname. 

STATUS 

Display  key  program  variables  states  or  values. 

TABULATE 

Tabulate  values  of  time-series  or  paired  data. 

TIME 

Set  a  time-window  for  time  series  data. 

SCO 

Resume  processing  subsequent  to  error. 

SAB 

Abort  (stop)  processing  subsequent  to  error. 

DSPtAY 

Initiate  the  DSPLAY  program. 
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**  Set  a  time  window 

Tl  T-2D  T 
*• 

Qot  tfio  *rsw*  dstiS 

QE  STG«RAWOB:/SCIOTO/HiGH3^AGE/01JUL198S/1  HOUR/OBS/ 

** 

**  Get  a  iook-up  tabie  for  the  station 

GE  TB=TABLEDB:/SCiOTO/HIGH3/STAGE-FLOW//AJSGS  TABLE  6/ 

** 

**  Compute  flows  using  the  lookup  table 

COM  FLW>=RrABLE(STG.TB) 

** 

**  Save  the  result 

PUT  FLW«MASTDB:/SCIOTO/HIGH3/FLOW/01JUL1 98^/1  HOUR/COMP/ 

** 

**  Contingency  checkpoint 

SCOmiNUE 

** 

**  Ciear  program  memories 
CL  ALL 


•*  Terminate 
Fi 

A  typical  command  sequence  retrieves  a  time  series  record,  computes  a  dependent  time 
series,  and  stores  the  result  An  example  input  sequence,  that  computes  flow  values  from 
stage  values  is  shown  in  Rgure  1. 


3^1  Data  Management 

3.2.1 .1  Program  Memory.  The  standard  memory  configuration  can 
accommodate  10  time  series  records  of  up  to  5300  values  each,  5  paired-data  records  each 
with  30  dependent  variables  with  a  maximum  of  5000  values  per  record,  and  15  scalars.  The 
DOS  version  configuration  is  limited  to  5  time-series  records  of  3180  values  each,  3  paired- 
data  records  eadi  with  5  dependent  variables  with  250  values  per  record,  and  15  scalars. 

Time  series  records  may  be  regular  or  irregutar  and  ail  data  times  are  stored.  Data 
quality  flag^  for  each  data  value  is  also  stored.  The  flags  indicate  four  levels:  1)  N  for  no 
flag;  2)  E  for  estimated;  3)  Q  for  questioned;  and  4)  M  for  undefined  or  missing.  Time  series 
records  may  be  retrieved  from  files  or  may  be  computed  in  DSSMATH.  However,  currently, 
data  quality  flags  may  not  be  read  or  stored  to  DSS. 
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Paired  data  recorda  include  rating  tables  Oe..  atage-flow)  arKi  polynofnial  coefficients. 
Paired-data  records  are  usually  retrieved  from  files,  although  at  the  present  time  only  one 
compute  function.  "MATE”,  creates  a  paired-data  record. 

Scalars  are  variables  that  are  results  ol  assignments  or  computations.  Storage  for 
scalars  is  provided  in  order  to  use  the  scalar  result  of  a  computation  in  a  subsequent 
computation. 

Each  record  or  variable  in  the  program  memory  is  assigned  an  aipharujmeric  label  as  it 
is  retrieved  or  computed.  The  label  Is  then  used  to  refer  to  that  time-series  in  further 
processing.  The  label  may  be  one  to  six  characters  long.  It  must  begbi  vvith  an  alphabetic 
character  but  may  contain  numeric  characters.  A  label  should  not  be  the  same  as  any  of  the 
commands  or  function  names. 

The  most  recently  used  time  wirKfow,  pathname,  and  DSS  fUe  name  are  kept  in  the 
program’s  memory.  A  maximum  of  five  DSS  flies  can  be  kept  open  on  a  continuous  basis. 
Once  a  sixth  DSS  file  is  requested,  the  first  DSS  file  opened  is  automaUcaily  closed  and  the 
new  DSS  file  requested  is  opened.  The  number  of  D^  files  that  can  be  open  on  a 
continuous  basis  is  iimited  to  three  fbr  the  640K  DOS  versioa 


3^1.2  Data  Retrieval  and  Storage.  The  commands  "GP  (get)  and  "Pir 
(put)  are  used,  respectively,  to  retrieve  arKl  store  time  series  and  paired  function  data.  The 
location  of  the  data  in  DSS  files  is  determined  by  both  a  file  name  and  a  pathname  (See 
Rgure  1).  Data  may  be  retrieved  and  stored  in  different  files  as  desired.  A  warning  message 
is  printed  if  the  retrievai  exceeds  any  storage  limits.  Options  are  provided  with  the  command 
"PIT  to  control  overwriting  of  previously  stored  data 

The  time  window  command  "IT  controls  the  span  of  the  data  retrieved  with  ”GP  and, 
therefore,  defines  the  period  of  time  for  processirrg.  The  time  window  must  be  defined  before 
any  processing  can  be  accomplished.  The  time  window  should  be  defined  carefully,  since 
arithmetic  computations  and  several  of  the  compute  functions  check  for  concurrence  of  the 
data 


When  regular-interval  time  series  data  are  retrieved,  the  time  window  is  temporarily 
adjusted,  if  necessary,  to  induda  data  belonging  to  exact  kitervai  boundaries.  Options  with 
the  retrieval  command  *GP  provide  retrievals  outside  the  range  of  the  time  window  to 
facilitate  use  of  irregular-interval  time  series  dafe. 


3.2.1. 3  Memory  Management  Functions.  The  status  command  fST) 
may  be  used  to  show  memory  states  and  contents,  including  ail  variable  label  names  and 
typ^  (time  series,  paired  function,  or  scalar),  the  current  time  window,  and  the  most  recently 
used  DSS  pathname  and  DSS  filename.  "ST”  options  will  also  show  the  pathnames,  units  arto 
types  of  time  series  and  paired-function  variables  and  values  assigned  to  scalars.  The 
command  "CL”  releases  memory  spaces  eitoer  globally  or  selectively,  by  label. 


3.2.2  Computations.  Computations  result  in  the  computation  of  values  of  a 
dependent  time  series,  paired  function,  or  scalar  from  a  simple  arithmetic  operation  or  a  more 
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complex  function  of  one  or  more  time  series,  paired  functions,  or  scalars.  If  the  dependent 
time  series  or  paired  function  variable  is  not  in  memory  it  is  created  with  undefined  pathname, 
type.  arKf  units.  However,  the  type  und  units  may  be  assigned  by  the  function  being  used. 

The  computation  of  a  depender.i  Ume  series  may  be  subjected  to  a  test  based  on  a 
logical  IF  corKfition  that  compares '  Jues  of  two  variables  or  a  variable  and  a  constant  The 
variables  may  be  concurrent  time  series  or  scalars,  if  the  IF  condition  is  not  satisfied  for  a 
particular  time,  then  the  corresporxfing  value  of  the  d^)endent  variable  is  unchanged.  If  the 
value  was  not  previously  defined.  It  wHi  rwnabi  urKiefined. 

In  general,  il  is  left  to  the  user  to  provide  appropriate  parameters  computations:  for 
example,  temperature  values  may  be  added  to  flows.  However,  arithmetic  computations  do 
require  concunent  time  series  variables,  and  DSS  data  types  or  units  are  checked  in  certain 
compute  functions.  Usually,  data  types  and  units  of  a  result  must  be  specified  explicitly  as  a 
separate  processing  step  with  the  *S0*  commarKi. 


3JL2.1  Arithmetic  Computations.  An  arithmetic  computation  of  a 
time-series  record  consists  of  addition  (+),  subtraction  (-),  multiplication  (*),  division  (/),  or 
exponentiation  (**)  by  a  constant,  previously  defined  scalar  or  time-series.  Tho  foliowing  are 
examples: 

COM  YaX4-2  add  2  to  time  aeries  *X* 

COM  YsX/A  "A*  may  be  a  acaiar  or  time-eeriea 

One  of  the  indeperxient  variables  in  an  arithmetic  computation  may  be  redefined  (ie.. 
the  dependent  variable  may  be  the  same  as  one  of  the  indeperxient  variables).  However,  if  a 
scalar  and  a  time  series  variable  are  the  independefft  variables,  the  dependent  variable  is  a 
time  series. 

When  arithmetic  computations  involve  two  hrKiependent  time-series,  the  dependent 
variable  will  consist  of  the  arithmetic  combination  of  the  concurrent  successive  values  of  the 
two  series.  "Concurrent  successive  values"  means  equal  numbers  of  values  with  the  same 
times.  The  "TSI"  or  "TS4"  compute  functions  are  useful  for  aligning  two  parallel  time  series 
vectors  in  time.  The  TSI"  function  generates  a  reguiar-intenrai  time  series  by  interpolation  at 
regular  intervals.  The  TS4"  function  generates  a  time  series  interpolated  to  the  times  of  a 
pattern  time  series. 

Two  time  series  involved  in  an  arithmetic  computation  may  have  different  units  and 
types.  Units  and  types  of  results  must  be  explidtiy  defined  by  the  user  with  the  "SO" 
command. 
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3.2.2.2  Functions.  Ck>mputations  may  also  Involve  special  functions  that 
typically  operate  upon  one  or  more  time-series  records  artd  result  In  a  scalar  or  new  time- 
series  record.  The  functions  are  summarized  in  Table  2  artd  described  an  appendix  B. 


Name 

Table  2  Compute  Functlona 

Oeacription 

ABS 

Absolute  function. 

ACC 

Running  accumulation. 

CMA 

Centered  moving  average  smoothing. 

CORR 

Compute  correlation  coefficients. 

COS 

Cosine  trigonometric  function. 

COUNTT 

Count  the  number  of  valid  and  missing  data. 

DDT 

Differences  per  unit  time. 

DECPAR 

Decaying  basin  wetness  parameter. 

DIFF 

Successive  differences. 

ESTUN 

Estimate  values  for  missing  data. 

ESTPPT 

Estimate  values  for  missing  predpitation. 

FMA 

Forward  moving  average. 

GENTSR 

Generate  a  regular  interval  time  series. 

INT 

Truncate  to  whole  numbers. 

LOG 

Natural  log  base  'e”. 

LOG10 

Log  base  10. 

MATE 

Generate  data  pairs  from  two  time  series  variables. 

MAX 

Maximum  value  in  a  time  series. 

MEAN 

Mean  value  in  a  time  series. 

MIN 

Minimum  value  in  a  time  series. 

MRG 

Merge  two  time  series. 

MUSK 

Muskingum  routing  function. 

NINT 

Round  to  nearest  whole  number. 

OLY  . 

Olympic  smoothing. 

PERCON 

Period  constants. 
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Name 

Table  2  (continued) 

Deecription 

POLY 

Polynomial  transformation. 

P0LY2 

Polynomial  transformation  with  integral. 

PULS 

Modified  Puis  or  Working  R&D  routing  function. 

QAC 

Flow  accumulator  gage  processor. 

RND 

Round  off. 

RTABLE 

Rating  tabie  interpoiation. 

RTABLR 

Reverse  rating  table  interpolation. 

RTABL2 

Two  variable  rating  tabie  interpoiation. 

SCRN1 

Scrcwtn  for  possible  erroneous  values  based  on  maximum  and 
minimum  range. 

SCRN2 

Screen  for  erroneous  data  vedu^  based  on  a  forward  moving  average 
maximum  value. 

SDEV 

Compute  the  standard  deviation  of  one  independent  variable. 

SHIFT 

Shfft  adjustment 

SIN 

Sine  trigonometric  function. 

SKEW 

Compute  the  skew  coefficient  of  one  independent  variable. 

SORT 

Square  root  function. 

SS 

Straddle  Stagger  routing  function. 

SSW 

Willmington  District  Straddle  Stagger  routing  function. 

TAN 

Tangent  trigonometric  function. 

TS1 

Interpolate  data  at  regular  time  intervals. 

TS2 

Period  averages  at  regular  intervals. 

TS3 

Period  minimums  and  maximums  at  regular  intervals. 

TS4 

Interpolated  data  at  irregular  intervals. 

TSCYCL 

Time  Series  Cydic  Analysis. 

TSHIFT 

Shift  time  series  in  time. 

TTSR- 

Transform  time  series  to  regular. 

TTSI 

Transform  time  series  to  irregular. 

1/X 

Inverse  function. 
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In  addition  to  the  functions  summarized  above,  DSSMATH  allows  the  sophisticated  user 
who  is  proficient  in  FORTRAN,  to  create  additional  compute  functions.  See  appendix  C  for 
Instructions  on  creating  additional  functions. 

Some  computations  result  in  a  scalar  vsyue.  Scalar  values  may  also  be  assigned  labels 
and  used  in  subsequent  processing,  but  cannot  be  retrieved  or  stored  in  OSS  files. 

3.2.3  Screening  and  Replacement 

The  screening  and  replacement  capabilities  of  DSSMATH  are  implemented  in  the 
compute  functions  ■SCRN1",  •SCRN2",  "ESTUN*,  and  •ESTPPT'.  Data  values  that  exceed  the 
spewed  limits  are  flagged  internally  to  facilitate  subsequent  corrections  with  estimation 
functions  or  graphical  editing. 

3.2.4  Plotting  and  Tabulation 

Plot  and  tabulation  capabilities  of  DSSMATHD  are  nearly  identical  with  those  of  the 
program  DSPLAY.  Time  series  and  paired  function  data  may  be  piotted  or  tabulated  with  the 
"PL”  and  TA*  commands,  respectively.  Commands  that  control  the  various  plotting  and 
tabulation  options,  such  as  grids  and  shading,  are  provided  in  DSSMATHD.  Refer  to  the 
DSPLAY  user’s  manual  for  additional  instructions  (HEC.  HEC-DSS  Users  Guide  and  Utility 
Program  Manualsl.  The  plotting  option  is  available  only  on  the  Harris  version.  The  tabulate 
option  is  available  on  all  versions  of  DSSMAim  and  is  no  longer  dependent  on  the  DSPLAY 
software. 

3.2.5  Miscellaneous  Commands 

The  DSS  file  catalog  command  "CA*  is  used  for  listing  the  pathnames  of  records  in  a 
DSS  file.  Data  variable  pathnames  and  descriptions  (eg.,  units)  can  be  specified  using  the 
*SP"  and  *SD"  commands,  respectively.  On-line  heip  is  available  through  the  "HE"  command. 
The  finish  command,  "FI",  terminates  processing. 


3.3  Contingency  Processing 

Certain  error  conditions  can  render  subsequent  processing  inappropriate.  For  example, 
if  data  cannot  be  retrieved,  it  cannot  be  processed.  These  errors,  which  simply  result  in 
warnings  in  the  interactive  environment,  may  be  processed  with  recovery  procedures  when 
input  are  taken  from  an  input  file  or  a  PREAD  macro  by  specifying  recovery  checkpoints  in 
the  input. 

An  input  iine  beginning  with  a  dollar  sign  ($)  is  an  error  recovery  checkpoint 
"$CONTINUP  indicates  where  to  resume  processing  when  an  error  is  encountered. 

"SABORT"  indicates  the  program  is  to  stop  if  an  error  has  been  encountered  and  no 
SCONTINUE  was  found,  if  no  checkpoints  are  found,  processing  stops. 
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Appendix  A 
DSSMATH  Commands 


Appendix  A  •  DSSMATH  Commands 
Command  Syntax 

Brackets  (Q)  and  ellipsis  (...)  symbols  are  notation  used  for  describing  commarKis  and  are  not 
part  of  the  actual  commarKi  syntax. 

UPPERCASE  items  in  uppercase  are  required  key  words  as  shown. 

lowercase  items  in  lowercase  are  to  be  supplied  by  the  user. 

[  ]  Items  within  brackets  are  optional. 

Items  immediately  preceding  ellipsis  may  be  repeated. 

Commas  and  blanks  are  used  to  separate  items  and  are  interchangeable  except  as  otherwise 
noted. 

Periods  are  required  when  options  are  specified.  No  blanks  should  appear  between  commands 
and  period  or  between  periods  and  options. 


Index  of  Commands 


Esas. 


**  Comment .  2 

CATALOG  Display  a  catalog  of  a  DBS  file  .  2 

CLEAR  Remove  a  variable  data  isdsei  from  memory .  3 

COMPUTC  Perform  a  computation .  3 

DIAG  Diagnostic  DSS  trace  .  4 

FINISH  Terminate  and  exit  the  program .  5 

GET  Retrieve  data  from  a  D^  file .  5 

HELP  List  commands .  6 

PLOT  Plot  variables  .  6 

PUT  Store  data  in  a  DSS  file .  6 

SD  Set  data  descriptions  .  7 

SHOW  Display  internal  data  vaHabie  information .  8 

SMISSING  Set  missing  value  indicators  .  8 

SP  Set  data  pathname .  9 

STATUS  Display  key  program  variables  states  or  values .  9 

TABULATE  Tabulate  values  of  time-series  or  paired  data .  10 

TIME  Set  a  time-window  for  time-series  data  .  10 

SCO  Resume  processing  subsequent  to  error .  10 

SAB  Abort  (stop)  processing  subsequent  to  error .  11 

DSPCAY  Initiate  the  DSPLAY  program .  11 
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Nam«: 


Comment 


Use:  **[parameters]  ... 

Description:  May  be  used  to  annotate  input  command  streams. 

Parameters: 

text  -  Any  text  up  to  1 30  characters  per  line. 

Example: 

**  Set  the  time  window 

Name:  CATALOG  Display  a  catalog  of  a  DSS  file 
Use:  CA  [.options]  [parameters] 

Description:  Catalog  (list)  the  records  (pathnames)  in  a  DSS  file.  The  catalog  is  listed  at 

the  terminal  one  screernfull  at  a  time.  The  numbers  associated  with  each 
record  may  be  used  in  other  commands  to  refer  to  the  record  in  lieu  of  a 
pathname. 

Options: 

None  (list  old  catalog  in  full  mode) 

N  (create  new  catalog) 

A  (create  abbreviated  catalog) 

Parameters: 

None  Catalog  last  opened  or  referenced  DSS  filename, 
filename:  Is  the  name  of  a  DSS  file. 

Example: 

CA.NA  MASDSS: 
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Name:  CLEAR 


Remove  a  variable  data  label  from  memory 


Uee:  CL  [parameters] 

Deecription:  dear  releases  memory  slots.  If  the  parameter  'ALL'  Is  presenL  all  slots  are 

released.  If  any  labels  are  specified,  then  only  the  memory  slots  represented 
by  the  labels  are  released. 

Parameters: 

None  No  action  Is  taken. 


ALL  All  variable  labels  are  cleared  and  all  data  is  initialized, 

label...  Will  dear  specified  variable  labels. 


Example: 

CLEAR  ALL 


Name:  COMPUTE  Perform  a  eomputaVonn 


Use:  CO  [parameters] 

Description:  DSS  data  header  information  are  checked  in  some  computations.  If  any  of 

the  independent  variable  values  are  flagged  as  missing  then  the  computed 
value  is  missing.  Data  types,  units  of  the  dependent  variable  are  generally 
undefined  but  may  be  set  by  some  functions. 

Parameters:  [IF(x1  operator  x2)l  v9sult»expression 

"resutr  is  a  labei  for  the  resultant  of  the  computation.  The  result  may  be  a  currently 
defined  variable  and  may  be  uexi  ?s  an  independent  variable  in  subsequent 
computations. 

"operator*  is  one  of  the  following  relationships: 

LT  less  than  LE  less  than  or  equal 

EQ  equal  NE  not  equal 

QT  greater  than  GE  greater  than  or  equal 


Appendix  A*3 
Commands 


*^1*  and  are  time  series,  scalars  w  constants.  Time  series  should  be  concurrent  with 
the  operand. 

'^expression”  is  either  a  simple  arithmetic  operation  or  a  function.  It  may  reference  one  or 
more  labels  for  variables  whk^  have  been  previously  computed  or  retrieved. 


Name:  DIAG  Diagnostic  DSS  trace 


Use:  Dl  [parameters] 

Description:  Toggles  diagnostic  trace  output  on  arKf  off.  Caution,  this  command  can 

generate  extensive  output. 

Parameters: 

None  Diagnostic  trace  is  off. 

ON  Turn  on  diagnostic  trace. 

OFF  Turn  off  diagnostic  trace. 

Example: 

01  ON 


Name:  DPATH  Display  a  selective  catalog  of  a  DSS  file 
Use:  OP  [.options]  [parameters] 

Description:  Displays  pathnames,  tags  and  reference  numbers  from  the  catalog  file  (see 

CATALOG  command).  The  DP  command  has  the  selective  display  capability 
which  provides  the  option  of  selecting  and  listing  only  certain  pathnames 
based  on  matching  pathname  parts. 

Options: 

None  (uses  last  referenced  DSS  file  catalog) 

N  (create  new  catalog) 

A  (create  abbreviated  catalog) 

Parameters: 

None  Displays  all  pathnames, 

filename:  Is  the  name  of  a  DSS  file. 

A=..B=..etc  Selective  catalog  based  on  pathname  parts. 

Example: 

DP  MASDSS:  B=RED  CREEK  C»FLOW 
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Name:  FINISH 


Tarmlnata  and  axit  tha  program 


Uaa:  FI 

Nama:  GET  Ratrlava  data  from  a  DSS  flla 

Uaa:  GE  [options]  [parameters] 

Description:  Retrieves  data  from  a  OSS  data  file.  The  data  may  be  time>series,  or  paired- 

function  data  ( eg.  rating  tables  or  polyriomlal  coefficients ). 

Options: 

P  Includes  the  value  just  prior  to  the  time  window. 

N  Includes  the  next  value  after  the  end  of  the  time  window. 

Parameters:  label=[filename:]pathname 

Taber  is  an  alphanumeric  identifier  for  the  data. 

"filename*  is  the  DSS  file  name  to  use  and  it  will  be  automatically  opened  if  necessary. 
If  omitted,  the  previously  opened  DSS  fie  is  used. 

"pathname*  may  be  an  explicit  pathname  or  a  pathname  catalog  number,  or  the  previously 
specified  pathname  may  be  modified  by  specifying  replacement  pathname 
parts.  If,  for  example,  the  previously  defined  pathname  were: 

/SCIOTO/CISG3/FLOW/01 FEB1 986/1 HOUR/OBS/ 

and  the  following  were  specified: 

B=HIGH3 

then  the  new  specification  would  become: 

/SCIOTO/HIGH3/FLOW/01 FEB1 986/1  HOUR/OBS/ 
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Name:  HELP  Ust  commands 


Use:  HE  [parameters] 

Description:  Displays  or>4rieclocumeritatioii  for  a  commarid.  If  rK)  command  is  given,  a  list 

of  commarKis  is  displayed.  In  order  to  get  a  list  of  the  functions,  enter  *HE 
FUNCTION". 

Parameters: 

"name*  Command  or  function  name  for  which  a  detailed  description  is  needed. 
FUNCTION  Will  display  a  listing  of  the  available  DSSMATH  functions. 


Name:  PLOT  Plot  variables 


Use:  PL  [parameters] 

Description:  Plots  the  data  represented  by  labeT  using  the  integrated  capabilities  of 

DSPLAY.  Up  to  seven  tedsels  may  be  plotted.  Both  time  series  and  paired 
function  data  may  be  plotted. 

Parameters: 

"laber  te  used  to  refer  to  the  data  v^in  DSSMATH.  Up  to  seven  different  labels  may 
be  specified. 


Name:  PUT  Store  data  In  a  DSS  file 
Use:  PU  [options]  [parameters] 


Description:  The  PUT  command  is  used  to  tidce  the  data  in  a  variable  label  and  store  it  in 

a  DSS  file  if  no  parameters  are  specified,  the  data  will  be  stored  in  the 
variables  pathname  if  one  has  been  given  or  in  the  last  pathname  defined. 

Options: 

A  All  data  in  the  time-series  identified  by  label  is  written  out  to  the  DSS  file  "filename" 
arKl  DSS  record  pathname,  replacing  any  existing  data. 

R  Replaces  existing  data  but  it  will  not  write  out  new  DSS  records  consisting  entirely 
of  missing  values. 

M  Same  as  option  R,  except  that  existing  data  will  not  be  replaced  by  a  missing 
value. 
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DeteuK  is  to  replaoe  only  missing  vakiss;  existing  non-missing  values  will  not  be 
replaced. 

Parameters:  label»[fllename:]pathname 

labeT  identifies  the  time  series  to  be  written 

"filename*  is  the  name  of  a  DSSfHe  to  receive  the  data.  If  it  is  not  specified,  then  the  last 
previously  referenced  "fiiename"  wUi  be  used.  The  DSS  file  will  be 
automatically  opened,  if  necessary. 

"pathname"  is  a  definition  of  the  pathname  to  use.  The  full  pathname  may  be  explicitly 
defined,  or  pathrwme  parts  may  be  i»ed  to  modify  the  last  previously  defined 
pathrrame  or  the  pathname  associated  with  the  label,  if  defirred. 


Name:  SD  Set  data  descriptions 

Use:  SO  [parameters] 

Description:  Data  descriptions  are  data  Items  contained  in  DSS  time  series  and  paired 

function  headers. 

Parameters:  label,parameter 

laber  is  the  DSSMATH  label  for  data. 

"parameter"  IrKiicate  vvhich  Rems  to  change  and  what  the  new  values  are.  "parameter" 
has  the  form  *Rem*s*value*.  where  "Hem*  is  TYPE  or  UNITS  for  time  series 
data  and  TYPE.  UNITS  or  LABELS  for  paired  function  data.  The 
designation  for  "Rem*  may  be  abbreviated,  such  as  U^  for  UNiTS^  . 
"value*  is  an  appropriate  entry  for  the  particular  "Rem*.  For  paired  function 
data,  enter  one  "value*  for  each  curve  in  a  series  separated  by  commas. 
For  example:  U=FEET,CFS. 
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Name:  SHOW 


Display  internal  data  variable  information 


Use:  SH  [parameters] 


Description:  Displays  selected  internal  information  about  the  data  variable  denoted  as  a 

parsuneter  to  the  SHOW  command.  Its  main  use  is  for  debugging.  However, 
this  command  is  the  only  way  to  see  the  value  of  a  scalar  variable.  For  time- 
series  K  shows  the  data's  array  position,  time  ( internal  representation),  time 
(external  representation),  value,  and  data  quality  flag.  For  paired  data,  it 
shows  the  DSS  header  information. 


Parameters:  label,parameter 


None  No  action  is  taken. 

label  Label  is  the  name  of  the  data  variable  that  will  be  displayed. 

SCALARS  Display  all  scalar  variables  and  their  values. 

Example: 


SH  STG 


Name:  SMISSING  Set  missing  value  Indicators 


Use:  SM  [parameters] 


Description:  Missing  value  indicators  are  used  to  define  the  numeric  values  with  which 

missing  data  is  indicated.  Up  to  10  numbers  can  be  specified.  If  no 
parameters  are  specified,  default  values  -901 .0  and  -902.0  are  used. 
Caution:  DSS  will  only  accept  the  detault  values  of  -901 .0  and  -902.0  as 
valid,  therefore  be  careftjl  about  writing  data  back  to  your  DSS  file  when  the 
missing  value  indicators  have  been  changed  from  the  default  settings. 

Parameters:  "parameter*  All  valid  numeric  values. 
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Name:  SP 


Sat  data  pathnama 


Uaa:  SP  [parameters] 


Daacription:  The  SP  commarid  is  used  to  define  or  set  the  default  pathname  for  a  variable 

l^)el  name. 

Paramatara:  label, parameter 

laber  is  the  DSSMATH  iabel  for  data. 

"parameter*  consists  of  either  a  pathname  or  a  patfmame  part  if  "parameter*  is  a  part  it 
has  the  form  "parr»*value*.  \where  "parf  is  A.  B.  C,  0.  E,  or  F,  and  "value*  is 
appropriate.  If  parts  are  specified,  the  resulting  pathname  for  the  data 
consists  of  the  cmrent  default  pathname  modified  by  replacement  of  the 
specified  parts. 


Name:  STATUS  Display  key  program  variables  states  or  values 


Use:  ST  [options]  [parameters] 

Description:  The  status  command  may  be  used  to  check  on  the  state  of  key  variables  and 

to  show  data  descriptions.  Those  functions  are  indicated  with  command 
options  and  parameters: 

Options: 

P  Display  pathname  for  variable. 

H  Display  data  descriptions  (DSS  header  contents). 

Parameters: 

VAR  Show  memory  allocations.  Icfoels,  and  associated  header  information. 

TIM  Show  the  current  time  window  setting. 

PA  Show  the  current  pathname. 

ALL  Aii  nf  the  above. 
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Nam#:  TABULATE 


Tabulata  valuaa  of  tima^rtas  or  palrad  data 


Uaa:  TA  [options]  [parametars] 

Daacrtption:  Tabulates  tha  data  represented  by  the  variable  label.  Up  to  seven  labels  may 

be  tabulated.  Both  time  series  or  paired  function  data  may  be  tabulated. 

Optiona: 

F  Send  the  output  to  the  file  specified  on  the  execution  line  by  the  use  of  the 
parameter  TAJBsfilename. 

Paramatars: 

Taber  Name  or  label  used  for  the  data  varii^. 

Nama:  TIME  Sat  a  tima^ndow  for  tlma-aarlaa  data 


Uaa:  Tl  [parameters] 

Description:  Starting  and  ending  times  and  dates  may  be  expressed  in  a  variety  of  ways, 

including  implicit  and  relative  times  and  dates: 

04MAR1 983  0700  Complete,  explicit  expression. 

D  0700  current  date,  time  expHcit 

D  current  date,  time  implicit  (2400) 

T  current  date  artd  time 

T-2H  two  hours  ago 

T-3D  three  days  ago 

T-15M  15  minutes  ago 

D-300  30  days  ago,  2400  hrs 

•0  +D  exparvds  the  exisfing  time  window  one  day  at  each  end 


Name:  SCO  Resume  processing  subsequent  to  error 


Use:  SCO 

Description:  Recognized  only  in  a  batch  job  or  when  a  PREAD  macro  is  being  used.  Used 

to  desigrtate  a  point  in  the  iri^  command  stream  to  resume  processing  when 
a  preceding  error  has  caused  processing  to  terminate. 
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Name:  $AB  Abort  (stop)  processing  subsequent  to  error 


Use:  SAB 

Description:  Recognized  only  in  a  batch  job  or  vwhen  a  PREAO  macro  is  being  used.  Used 

to  specify  termination  of  ttie  job  if  an  error  is  encountered  while  processing 
according  to  a  batch  input  stream. 


Name:  DSPLAY  initiate  the  DSPLAY  program 


Use:  DSPLAY 

Description:  The  DSPLAY  command  will  initiate  the  graphics  program  DSPLAY.  The  user 

can  either  plot  pathnames  from  the  latest  DSS  file  opened  in  DSSMATH  or  the 
user  can  use  the  OPEN  command  to  open  another  DSS  file.  The  majority  of 
the  regular  DSPLAY  commands  are  available.  To  exit  the  DSPLAY  program, 
the  user  must  enter  the  FINISH  commarKf.  (This  command  is  only 
available  on  the  Harris  computer  under  program  DSSMATHD) 
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Appendix  B  -  Compute  Functlone 


Index  of  Functlone 

Page 

ABS  Absolute  function .  3 

ACC  Running  accumulation  .  3 

CORR  Compute  correlation  coefficients .  3 

COS  Cosine  trigonometric  function .  4 

COUhiT  Count  the  number  of  valid  and  missing  data  values  .  4 

CMA  Centered  moving  average  smoothing .  5 

DDT  Differences  per  unit  time .  6 

DECPAR  Decaying  basin  wetness  parameter .  6 

DIFF  Successive  differences .  6 

ESTUN  Estimate  values  for  missing  data .  7 

ESTPPT  Estimate  values  for  missing  precipitation  data .  7 

FMA  Forward  moving  average  .  7 

GENTTSR  Generate  a  regular  interval  time  series  .  8 

ihn*  Truncate  to  whole  numbers  .  8 

LOG  Natural  log  base  "e*  .  8 

LOG10  Log  base  10 .  8 

MATE  Generate  data  pairs  from  two  time-series  .  9 

MAX  Maximum  value  in  a  time  series  .  9 

MEAN  Mean  value  in  a  time  series .  9 

MIN  Minimum  value  in  a  time  series .  9 

MRG  Merge  two  time  series .  10 

MUSK  Muskingum  hydrologic  routing .  10 

NINT  Round  to  nearest  whole  number .  10 

OLY  Olympic  smoothing .  11 

PERCON  Period  constants .  11 

POLY  Polynomial  transformation .  12 

POLY2  Poi^mial  transformation  with  integral .  12 

PULS  Modified  Puls  or  Working  RU)  routing  function .  12 

QAC  Flow  accumulator  gage  processor .  13 

RND  Roundoff .  13 

RTABLE  Rating  table  interpolation  .  14 

RTABLR  Reverse  rating  table  inteipoiation  .  14 

RTABL2  Two^riable  rating  table  interpolation .  14 

SCRN1  Screen  for  possible  enoneous  values  based  on  maximum/minimum 

range .  15 

SCRN2  Screen  for  possible  erroneous  values .  15 

SDEV  Compute  the  standard  deviation  of  one  independent  variable .  15 

SHIFT  Shift  adjustment .  16 

SIN  '  Sine  trigonometric  function .  16 

SORT  Square  root  function .  16 
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Functions 


ss 

ssw 

TAN 

TS1 

TS2 

TS3 

TS4 

TSCYCL 

TSHIFT 

TTSR 

TTSI 
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Straddle  Stagger  routing  function .  17 

WiUmington  District  Straddle  Stagger  routing  function .  17 

Tangent  trigoiiornetric  furtction .  17 

Interpolated  data  at  regular  intervals .  18 

Period  averages  at  regular  intervals  .  18 

Period  mins  or  maxs  at  regular  intenrals .  18 

Interpolated  data  at  irregular  intervals .  .  19 

Time  series  cyclic  analysis  .  19 

Shift  time  series  in  time .  20 

Transform  time  series  to  regular .  20 

Transform  time  series  to  irregular .  21 


Name: 


ABS 


Absolute  function 


Use: 

Description: 

Name: 

Use: 

Descripliofi! 

Name: 

Use: 

Description: 


CO  TY=ABS(TX) 

Compute  the  absolute  value  of  TX  and  store  in  TY.  TX  and  TY  can  be  the 
same  variable.  TX  can  be  soalar,  ttone  series  or  paired  data.  IfavalueofTX 
is  undefined  resulting  values  of  TY  are  undefined.  If  a  concurrent  IF 
condition  is  rrat  satisfied,  TY  is  unchanged. 


ACC  Running  accumulation 


COTY=ACC(TX) 

Compute  a  running  accumuiation  of  the  values  in  TX  and  store  the  result 
in  TY.  If  a  TX  value  is  undefined  or  a  concurrent  IF  condition  is  not  satisfied, 
the  value  of  TX  is  rK>t  added  to  the  accumulation  and  the  corresponding  TY 
remains  the  same  as  the  previous  TY.  Units  of  TY  are  the  same  as  those  of 
TX  If  TX  is  typed  as  PERAVER  or  PER-CUM,  TY  will  be  typed  INST-VAL  or 
INST-CUM,  respectively. 


CORR  Compute  correlation  coeffieients 


CO  TS«CORR(TX1  .TX2.INDEX) 

Compute  the  correlation  coefficients,  determination  coefficients,  and 
standard  errors  of  regression  between  the  two  variables  TX1  and  TX2.  TX1 
and  TX2  must  be  of  the  same  type  (either  uniform  time  series  or  irregular 
time  series)  and  contain  the  same  number  or  rows  or  ordinates.  Variabie 
INDEX  is  used  to  set  the  scalar  variable  TS  as  follows: 

Index  Description 

1  Number  of  valid  pairs  for  correlation. 

2  Regression  constant 

3  Regression  coefficient 

4  Determination  coefficient 

5  Standard  error  of  regression. 

6  Determination  coefficient  adjusted  for  degrees  of  freedom. 

7  Standard  error  adjusted  for  degrees  of  freedom. 
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Name: 


COS 


Coalna  trigonometric  function 


Use: 

Description: 

Name: 

Use: 

Description: 


CO  TY«COS(7X) 

Compute  the  coeine  of  TX  and  store  in  TY.  TX  and  TY  can  be  the  same 
variable.  TX  can  be  scalar,  time  series  or  paired  data.  If  a  value  of  TX  is 
undefined  resulting  values  of  TY  are  urKiefined.  if  a  cor^rrent  IF  condition 
is  not  satisfied,  TY  is  unchanged. 


COUNT  Count  the  number  of  valid  and  missing  data  values 


CO  TS=COUNT(TXTYPE) 

This  function  is  used  to  count  the  valid  and  missing  values  in  time  series 
variable  TX  and  stores  it  in  scalar  variable  TS.  Possible  values  for  TYPE  are: 
VAUD*  and  ‘MISSING*. 
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Name: 


CMA  Cantered  moving  average  emoothing 


Uae: 

Deacriptlon: 


LEVEL1 

LEVEL2 

LEVEL3 

LEVEL4 

Note: 


CO  TY-CMA(TX.NP[.CLEVEL]) 

Compute  a  centered,  moving  average  of  NP  values  in  TX  arKl  store  in  TY. 
NP  must  be  odd  and  greater  than  2.  The  default  CLEVEL  is  'LEVELa*  and 
NP/2  values  at  the  beginning  and  end  of  TY  will  be  urKiefined.  If  a  value  in 
TX  is  undefined  or  a  concurrent  IF  corKlition  is  not  satisfied,  then  the 
resulting  TY  values  are  the  averages  of  orte  less  value.  TX  must  be  a 
regular-intervai  time-series  and  TY  and  TX  cannot  be  the  same.  Usefcjifbr 
filtering  instantaneous  data  containing  high-frequency  variations.  Units  and 
type  are  unchanged.  The  following  levels  are  available: 

Only  valid  values  will  be  averaged.  Valid  values  that  do  not  have  all  valid 
values  within  the  averaging  period  number  will  be  missing.  Values  at  the 
start  and  end  of  the  data  that  do  not  have  enough  valid  values  to  average 
will  be  assigned  missing  values. 

Only  valid  values  will  be  averaged.  Valid  values  that  do  not  have  all  valid 
values  within  the  averaging  period  number  will  be  missing.  Valid  values  at 
the  start  and  end  of  the  data  that  do  not  have  enough  valid  values  to 
average  will  be  averaged  over  a  reduced  number  of  values. 

(Defoult  setting)  All  values  will  be  averaged  based  on  valid  values  within  the 
averaging  period  number.  If  all  values  are  missing  for  the  period  number, 
the  value  will  be  missing.  Values  at  the  start  and  end  of  the  data  that  do  not 
have  enough  valid  values  to  average  will  be  assigned  missing  values. 

All  values  will  be  averaged  based  on  valid  values  within  the  averaging  period 
number.  If  all  values  are  missing  for  the  period  number,  the  value  will  be 
missing.  Values  at  tiie  start  and  end  of  the  data  that  do  not  have  valid 
values  to  average  will  be  averaged  based  on  a  reduced  number  of  values. 

Questionable  and  estimates  flagged  values  are  used  in  the  computations. 
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Name: 


DDT 


Differancaa  par  unit  time 


Use: 

Description: 


Name: 

Use: 

Description: 


Name: 

Use: 

Description: 


COTY=DDTfTX) 

Computes  the  successive  differences  in  TX  per  time  period: 
7Y(t)=aX(t)-TX(t.1))/DT 

where  OT  is  the  time  difference  in  days  between  t  and  M.  if  a  value  of  TX 
is  undefined  TY  is  undefhed.  if  a  concurrent  IF  condition  is  not  satisfied,  TY 
isunchanged.  TY  and  TX  cannot  be  the  same  variable.  TY  is  assigned  the 
type  'PER-AVER.*  The  units  of  TY  are  urKiefined.  An  example  of  the  use  of 
DOT  is  the  computation  of  reservoir  inflow  from  outflow  and  the  change  in 
storage,  where  the  change  in  storage  is  transformed  to  average  flow  volume 
per  day  by  the  function.  The  conversion  factor  fc>r  storage  to  flow  is  0.5041 6 
when  the  storage  change  is  in  ao4t  day. 


DECPAR  Decaying  basin  wetness  parameter 


CO  TY=DECPAR(TY.T2.R) 

Compute  a  time-series  of  parameters  TY  as  a  furrcQon  of  TZ  and  R: 
TY(t)=:RTY(M)+TZ(t) 

where  R  is  the  decay  rate  and  TZ  is  precipitation.  R  is  less  than  1.  The 
function  extends  TY:  the  first  value  in  the  series  TY  is  used  as  the  starting 
value,  and  any  other  TY  are  computed  in  the  sequence.  If  the  first  value  in 
TY  is  undefined,  it  is  assumed  to  be  zero.  The  first  TZ  value  is  ignored.  TY 
and  TZ  must  be  regular-interval  time  series  with  identical  times.  R  should  be 
appropriate  for  the  intenral.  The  type  arfo  units  of  TY  are  unchanged.  An 
IF  coTKlition  has  no  effect 


DIFF  Successive  differences 


COTY=DIFF(TX) 

Compute  the  successive  differences  of  TX  and  store  in  TY.  TX  and  TY 
cannot  be  the  same  variable.  TX  must  be  of  type  INST-VAL  or  INST-CUM. 
If  a  value  of  TX  is  undefined,  resulting  values  of  TY  are  undefined.  If  a 
concurrent  IF  condition  is  not  satisfied,  TY  is  unchanged.  The  type  and 
units  of  TY  are  undefined. 
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Nam*: 


Ettimato  valuaa  for  mlaaing  data 


Uae: 

Deaerlption: 

Name: 

Uae: 

Deacription: 


Name: 

Uaa: 

Deacription: 


CO  TY«ESTUN(TXNMAX.QFLAQ) 

Unearly  Interpoiata  astbnates  for  values  in  TX  with  flags  equal  or  lesser  in 
quality  to  QFLAQ  and  place  the  residts  in  TY.  Do  not  estimate  more  that 
NMAX  continuous  missing  values.  Possft>le  flags,  hi  order  of  decreasing 
quality,  are:  Q  «  questtOTAbie  and  M  »  missing.  An  IF  condition  has  no 
effect  TX  must  be  a  time-series,  arKf  TX  arKf  TY  may  be  the  same  variable. 
Type  and  units  of  TY  are  the  same  as  TX 


ESTPPT  Eatlmata  valuaa  for  missing  praelpltation  data 


CO  TY=ESTPPT(TXNMAXQFLAG) 

Unearly  interpolate  estimates  for  cumulative  precipitation  values  in  TX  with 
quality  flags  equal  to  or  lesser  in  quality  to  QFLAG  and  place  the  results  in 
TY.  Possible  flags,  in  order  of  decreasing  quality,  are:  Q  »  questionable 
and  M  •  missing.  If  the  values  bracketing  the  missing  period  are  increasing 
with  time,  do  not  estimate  more  that  NMAX  oonthmous  missing  values,  if 
the  values  bracketing  the  missing  period  are  equal,  then  estimate  any 
number  of  missing  values,  if  the  values  bracketing  the  missing  period  are 
decreasing  with  time,  do  not  estimate  any  missing  values.  TX  must  be  a 
data  type  of  INST-CUM.  An  IF  corKlition  has  no  effoct 


FMA  Forward  moving  average 


C0TY=FMA(TXNP) 

Compute  a  moving  average  of  the  last  NP  values  in  TX  and  store  in  lY.  NP 
must  greater  than  2.  NP  values  at  the  beginning  of  TY  will  be  missing.  Ifa 
value  in  TX  is  missing,  the  value  is  not  used  for  computing  TY,  and  the 
average  is  over  one  less  value.  At  least  2  values  of  TX  must  be  defined, 
else  TY  is  missing.  Useful  for  computing  flow  durations.  Also,  may  be  used 
to  determine  parameters  for  use  in  the  SCR2  screening  function.  Units  and 
type  are  unchanged. 
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Nam*: 


QENTSR  Genarate  a  ragular  Intarval  lima  aariaa 


Uaa: 

Daacription: 

Nama: 

Uaa: 

Daacription: 

Nama: 

Uaa: 

Description: 

Name: 

Use: 

Description: 


CO  7Y=GENTSR(DT.TOFF.YO,QFLAQ) 

Generate  a  new,  regular-intofval  time  series  TY  with  a  time  krtervaJ  ot  DT,  a 
time  Intenral  offset  of  TOFF,  a  constant  value  and  an  values  internally 
flagged  with  QFLAG.  TOFF  is  tbne  from  the  beginning  of  the  standard 
interval  to  the  actual  time  of  the  data.  DT  and  TOFF  are  expressed  as  units 
of  time  as  rmT,  where  T  Is  M  for  minutes,  H  for  hours,  and  D  for  days. 
Possible  flags  are:  N  «  none,  M  «  missing.  E  «  estimated,  and  Q  > 
questionable.  lfY0is-801.,thentheflagisautomaticaliy M.  Unitsandtype 
must  be  set  indeperKlentiy  (use  SDD).  An  IF  condition  has  no  effect 

INT  Truncate  to  whole  numbara 


CO  TY=INT(TX) 

Truncate  to  a  whole  number  TX  and  store  In  TY.  TX  arKJ  TY  can  be  the 
samevariable.  TX  can  be  scalar,  time  series  or  paired  data.  IfavaiueofTX 
is  undefined  resulting  values  of  TY  are  undefined.  If  a  concurrent  IF 
coTKiition  is  not  satisfied.  TY  is  unchanged. 


LOG  Natural  log  base  *a* 


CO  TY=LOGCTX) 

Compute  the  natural  log  base  e  ofTXand  store  in  TY.  TXandTYcanbe 
the  same  variable.  TX  can  be  scalar,  time  series  or  paired  data.  If  a  value 
of  TX  is  UTKiefined  resulting  values  of  TY  are  undefined.  If  a  concurrent  IF 
condition  is  not  satisfled,  TY  is  unchanged.  If  a  value  in  TX  is  0  or 
negative,  the  value  in  TY  will  be  set  to  missing. 


LOGIC  Log  base  10 


COTY=LOG10(TX) 

Compute  the  log  base  10  of  TX  and  store  in  TY.  TX  and  TY  can  be  the 
samevariable.  TX  can  be  scalar,  time  series  or  paired  data.  IfavaiueofTX 
is  undefined  resulting  values  of  TY  are  undefined.  If  a  concurrent  IF 
condition  is  not  satisfied,  TY  is  unchanged.  If  a  value  in  TX  is  0  or 
negative,  the  value  in  TY  will  be  set  to  missing. 
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N&.*ne: 


MATE 


G«n«rato  data  paira  from  two  tlm»«eriaa 


Use: 

Doacription: 

Name: 

Use: 

Description: 

Name: 

Use: 

Description: 

Name: 

Use: 

Description: 


CO  PF=MATECTXTY.SORr/NOSOfn} 

Derive  a  paired-function  PF  by  pairing  values  in  the  time  series  TX  and  TY. 
TX  artd  TY  values  must  have  identical  times.  (Functions  TS1  arKi  TS4  may 
be  useful  in  corijunction  with  MATE).  The  paired  data  are  sorted  into 
ascerKlingorderofTXifSORTisspedlied.  The  urUts  of  PF  are  respective 
the  same  as  those  of  TX  and  TY.  The  types  of  PF  are  undefirwd.  An  iF 
coTKlition  has  no  effect 


MAX  Maximum  tmlua  In  a  time  series 

CO  SY=MAX(TX) 

Rnd  the  maximum  value  in  time-series  TX  and  place  it  in  scalar  SY.  Ignore 
missing  values  or  values  concurrent  with  an  unsatisfied  IF  corKfition. 

MEAN  Mean  value  In  a  time  series 

CO  SY=MEANCrX) 

Compute  the  mean  value  in  time-series  TX  and  place  the  result  in  scalar  SY. 
Ignore  missing  values  or  values  concurrent  with  an  unsatisfied  IF  condition. 

MIN  Minimum  value  In  a  time  series 


CO  SYsMIN(TX) 

Rnd  the  minimum  value  in  time-series  TX  and  place  the  result  in  scalar  SY. 
Ignore  missing  values  or  values  concurrent  with  an  unsatisfied  IF  condition. 
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Nairn: 


MRQ 


Marga  two  tlma  aariaa 


Uaa: 

Daacriptlon: 

Name: 

Uaa: 

Deacription: 

Name: 

Use: 

Description: 


CO  TY«MRG(TX.7Z.QFLAG) 

Merge  TX  with  7Z  TY  indudes  alt  values  In  TX  and  TZ,  except  where  TX 
arKiTZ  occur  at  the  same  time.  In  that  case,  the  value  for  TX  is  used  unless 
it  is  flagged  with  a  qualit/  equal  to  or  less  than  QFLAG  and  TZ  is  flagged 
with  a  quality  greater  than  QRAG.  Possible  flags,  In  order  of  decreasing 
quality,  are:  N  «  no  ftag,  E  «  estimated,  Q  «  questionable  and  M  « 
missing.  An  IF  oorKlition  has  no  effect  The  type  and  imits  of  TY  are 
undefined. 


MUSK  Muskingum  hydrologic  routing 


CO  TY»MUSK(TXNR,K,X) 

Route  the  uniform  time  series  variable  TX  by  the  Muskingum  hydrologic 
routing  method  and  store  it  in  variable  TY.  The  *IP  compute  options  has  no 
effect  on  this  function.  The  1C  parameter  is  the  Muskingum  V  in  hours,  1C 
parameter  is  the  Muskingum  V  (range  between  0  and  .5),  and  I^R*  is  the 
number  of  routing  subreaches. 


NINT  Round  to  nearest  whole  number 


CO  TY=NINT(TX) 

Round  to  nearest  whole  number  TX  and  storeinTY.  TXandTYcanbe  the 
same  variable.  TX  can  be  scalar,  time  series  or  paired  data.  If  a  value  of  TX 
is  undefined  resulting  values  of  TY  are  undefined.  If  a  concurrent  IF 
condition  is  not  satisfied,  TY  is  unchanged. 
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Olympic  tmoolhlng 


Name: 

Use: 

Description: 

LEVEL1 

LEVEL2 

LEVELS 

LEVEL4 

Note: 

Name: 

Use: 

Description: 


OLY 


CO  TY»OLYrX.NP[CLEVELI) 

Compute  a  smoothed  time-series  from  TX  using  the  Olympic  srrKXJthing 
scheme:  same  as  centered,  moving  average  except  the  minimum  and 
maximum  values  In  the  span  NP  are  igrKHed.  Place  the  resirit  in  TY.  Units 
and  type  are  imchanged.  The  following  levels  are  available: 

Only  valid  values  will  be  averaged.  Valid  values  that  do  not  have  all  valid 
values  within  the  averagirrg  period  number  will  be  missing.  Values  at  the 
start  and  end  of  the  data  that  do  not  have  enough  valid  values  to  avmage 
will  be  assigned  missing  vidues. 

Only  valid  values  will  be  averaged.  Valid  values  that  do  not  have  all  valid 
values  witNn  the  averaging  period  number  will  be  missing.  Valid  values  at 
the  start  arKf  end  of  the  data  that  do  not  have  enough  valid  values  to 
average  adll  be  averaged  o\m  a  reduced  number  of  values. 

(Oetault  setting)  AH  values  will  be  averaged  based  on  valid  values  within  the 
averaging  period  number,  if  all  values  are  missing  for  the  period  number, 
the  value  will  be  missing.  Values  at  the  start  and  end  of  the  data  that  do  not 
have  enough  valid  values  to  average  will  be  assigned  missing  values. 

All  values  will  be  averaged  based  on  valid  values  within  the  averaging  period 
number.  K  all  values  are  missing  for  the  period  number,  the  value  will  be 
missing.  Values  at  the  start  and  end  of  the  data  that  do  not  have  valid 
values  to  average  will  be  averaged  based  on  a  reduced  number  of  values. 

Questlonabie  and  estimates  flagged  values  are  used  in  the  computations. 


PERCON  Period  constants 


CO  TY=PERCON(rS,TX) 

Generates  a  series  TY  at  times  concurrent  with  TX  and  with  values  equal  to 
the  previous  chronological  TS  value.  If  no  previous  TS  is  present,  then  TY 
values  are  undefined.  TS  and  TX  may  be  irregular-  or  regular-intenml  time 
series.  An  IF  condition  has  rK>  effect 
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Nam«: 


POLY 


Polynoinlal  transformation 


Uaa: 

Daacrlption: 

Name: 

Use: 

Description: 

Name: 

Use: 

Description: 


COTY-POLY(nC.TP) 

C>>n^)ute  a  polyrKNTiial  transformation  of  TX  using  the  polynomial  coefficients 
TP.  Store  the  result  In  TY.  If  a  TX  value  is  missing  TY  is  undefined.  If  a 
concurrent  IF  condition  is  not  satisfied,  lY  is  unchanged.  Units  and  type  of 
TY  are  defined  by  TP.  TP  can  be  created  with  the  utility  DSSPO. 


POLY2  Poiynomiai  transformation  with  integrai 


COTY-*POLY2(TXTP) 

Compute  a  poiynomiai  transformation  of  TX  using  the  integral  of  the 
poiynomiai  defined  by  coefficients  TP.  Store  the  result  in  TY.  If  a  TX  value 
is  missing  TY  is  urKfefined.  if  a  concurrent  IF  condition  is  not  satisfied,  TY 
is  unchanged.  Units  arKi  type  are  uncharrged.  TP  can  be  created  with  the 
utility  OSSPO. 

PULS  Modifiad  Pula  or  Working  R&D  routing  function 


CO  TY«PULSCTX,PFX.NR,STOR1 ,001) 

Route  the  uniform  time  series  variable  TX  by  the  Modified  Puls  or  Working 
R&O  hydrologic  routing  method  and  store  it  in  variable  TY.  The  variable  PF 
must  have  been  previously  defined  in  a  GET  command.  It  references  a 
storage-discharge  paired  function  relationship  for  the  reach.  Storage  must 
be  the  first  variable,  discharge  the  second  variable,  and  each  variable  must 
repeat  only  once.  Variable  X  is  the  wedge  coefficient  (Muskingum  X)  kM’  use 
in  working  R&D.  Use  0.0  value  to  route  by  Modified  Puls  method.  Variable 
NR  is  the  number  of  routing  reaches.  STOR1  and  001  are  initial  storage 
and  flow  respectively.  Use  a  value  of  -1  for  both  STOR1  and  Q01  in  order 
to  use  the  first  flow  value  in  TX  and  interpolate  a  corresprxiding  storage 
value  from  PF.  The  IF  compute  command  has  no  effect  on  this  function. 


ApperKlix  B-12 
Functions 


Name:  QAC  Flow  accumulator  gaga  procosaor 


Uaa:  CO  TY«QAC(TX.TC) 

Daacription:  Compute  period-average  flows  from  a  flow  accumulator  type  gage: 

TY(t)  -  (TXCt^-DCtt-lM/frCW-TCCM)) 

TX  arKl  TC  are  respectively,  time  series  of  the  accumulated  flow  aiKl  the 
count  TX  and  TC  values  must  occur  at  the  same  times.  If  corresponding 
values  for  TX  atKJ  TC  decrease  from  the  previous  period,  thwi  the  accumula¬ 
tion  is  assumed  to  have  reset  to  zero  at  the  beginning  of  the  interval.  An  IF 
condition  has  rto  effect  TY  Is  assigned  the  type  'PER-AVER.' 


Name:  RND  Round  off 


Use:  COTY»RNO(TXNDIG,IPLAC) 

Description:  Round  off  values  to  NOIG  or  IPLAC,  whichever  controls.  NDIG  is  the 

number  of  digits  to  round  to  and  can  range  from  1  to  8.  IPLAC  is  a 
magnitude  of  1 0  to  which  to  round  to:  for  exampie,  -1  specifies  rourufing  to 
one-tenth  (0.1).  The  number  of  digits  shown  Is  never  less  than  1,  however. 
The  fbliowing  example  illustrates  the  effects  on  rounding  if  NDIG^S  and 
IPLACE=0  (ie.,  round  to  ones  place): 

rounds  to 

1445.1  •  1450. 

144.51  "  145. 

14.451  "  14. 

1.4451  •  1.0 

If  TX  is  undefined,  TY  is  undefined.  H  a  concurrent  IF  condition  is 
unsatisfied,  TY  is  unchanged. 
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Functions 


Name: 


RTABLE  Rating  table  Interpolation 


Uae: 

Deaerlptlon: 

Name: 

Uae: 

Description: 

Name: 

Use: 

Description: 


CO  7Y«RTABLE(D(,TB) 

Interpolate  values  fbr  TX  using  table  TB  and  store  the  resuK  In  TY.  TB  must 
be  created  using  the  program  DSSPD  (use  /R  option)  with  spedfic 
information  in  the  header:  type  of  interpola^.  ofteet.  shift,  and  datum,  if 
the  type  of  interpolation  is  LOGLOG,  table  x  values  are  adjusted  by 
subtracting  the  ofteet  The  shift  is  added  to  and  the  datum  subtracted  from 
all  incoming  TX  values.  The  header  information  in  TB  also  defines  the  units 
ofTY. 


RTABLR  Reverse  rating  table  Interpolation 


COTY*RTABLR(TX.TB) 

Interpolate  values  for  TX  using  the  reverse  of  table  TB  and  store  the  result 
in  TY.  TB  must  be  created  using  the  program  DSSPD  (option  /R)  with 
specific  information  in  the  header:  type  of  interpolation,  offset,  shift,  and 
datum,  if  the  type  of  inteipoiation  is  LOGLCXa,  table  x  values  are  adjusted 
by  subtracting  the  offset  The  shift  is  subtracted  from  and  the  datum  added 
to  ail  resulting  TX  values.  The  header  information  in  TB  also  defines  the 
units  of  TY. 

RTABL2  TWo<variabls  rating  table  Interpolation 


CO  TY=RTABL2(7X,TZ.TB) 

TY  is  a  function  of  two  independent  variables  TX  and  TZ.  The  functional 
relationship  is  specified  by  the  table  TB,  which  consists  of  sets  of  TX/TY 
pairs,  each  set  corresponding  to  a  TZ  value.  TB  is  created  with  the  program 
DSSPD  (option  /R)  with  values  fOr  TZ  specified  as  labels  for  the  sets  of 
TX/TY  pairs.  RTABL2  interpolates  linearly  in  table  TB.  No  extrapolation  is 
done:  if  the  TX  or  TY  value  are  outside  the  range  bounded  by  TB,  TY  is  set 
to  missing.  TX  and  TZ  must  be  concurrent  time  series.  The  header  informa¬ 
tion  in  TB  also  defines  the  unite  of  TY. 
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Name: 

Use: 

Description: 

Name: 

Use: 

Description: 

Name: 

Use: 

Description: 


SCRN1  Screen  for  possible  erroneous  values  based  on 
maximum/minimum  range 


CO  TY«SCRN1  (TXXMIN.XMAX.MAXDEL.QFLAG) 

Flag  any  value  bi  TX  that  tails  outside  the  range  XMIN  -  XMAX  or  exceeds 
the  maximum  charrge  MAXOEL  from  the  previous  value.  The  maximum 
change  comparison  is  done  ot^  when  the  consecutive  values  are  not 
flagged.  Possible  values  for  QFLAQ  are:  M  «  missing  data  or  Q  « 
questionable.  The  result  is  placed  in  TY. 


SCRN2  Screen  fw  possible  erroneous  values  based  on 
forward  moving  average  maximum 


CO  TY«SCRN2(TXNPTS.MAXDEL.QFLAQ) 

Flag  any  value  in  TX  that  exceeds  the  maximum  change  MAXOEL  from  the 
forward  moving  average  of  NPTS  ending  at  the  previous  value.  Missing 
values  in  TX  are  not  counted  In  the  moving  average  and  the  divisor  of  the 
average  is  less  orw  for  every  missing  value.  Values  which  fail  the  screen  are 
not  counted  either.  At  le^  2  values  must  be  defined  else  the  moving 
average  is  undefined  and  the  screen  passes  the  relevant  TX  value.  Possible 
flags  are:  M  missing  data  or  Q  «  questionable.  The  result  is  placed  in 
TY.  SCRN2  is  useful  for  detecting  and  removing  spikes.  The  FMA  function 
may  be  useful  for  determinlr^  appropriate  NPTS  and  MAXOEL  parameters. 


SDEV  Compute  foe  standard  deviation  of  one  independent 
variable. 


CO  TS=SOEV(TX) 

This  function  is  used  to  compute  the  standard  deviation  of  one  independent 
variable.  TX  can  be  either  regular  or  irregular  time  series.  TS  is  a  scalar 
and  is  set  to  the  computed  standard  deviation.  This  function  requires  a 
minimum  of  three  valid  values  to  compute. 
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Name: 


SHIFT 


Shift  adluatment 


Use: 

Description: 

Name: 

Use: 

Description: 

Name: 

Use: 

Description: 

Name: 

Use: 

DescriptlOii: 


CO  7Y«SHIFr(TB.TX) 

Generate  a  time  eerlee  cl  shift  adjustments  TY  with  times  at  TX  and  actuai 
shifts  TS,  possibly  at  other  limes.  TX  arKi  TY  onnot  be  the  same  variable. 
TS  values  are  interpolated  at  TX  times.  The  interpolation  is  linear  between 
TS  values,  ttcoept  when  TX  time  is  greater  then  last  TS  time.  Then,  the  last 
value  of  TS  is  held  constant  for  the  remaining  TX  times.  If  no  previous  TS 
value  brackets  TY,  then  TY  is  set  to  zero.  An  IF  condition  has  no  effect 

SIN  Sine  trigonometric  function 


CO  TY*SIN(TX) 

Compute  the  sine  of  TX  and  store  in  TY.  TX  and  TY  can  be  the  same 
variable.  TX  can  be  scalar,  time  series  or  paired  data.  If  a  value  of  TX  is 
undefined  resulting  values  of  TY  are  undefined.  If  a  concurrent  IF  condition 
te  not  satisfied,  TY  is  unchanged. 


SKEW  Compute  the  skew  coefficient  of  one  independent 
variable 


CO  TS=SKEW(TX) 

This  foncfton  is  used  to  ccNupute  the  skew  coefficient  of  one  independent 
variable.  TX  can  be  either  regular  or  irregular  time  series.  TS  is  a  scalar 
and  is  set  to  the  computed  skew  coeffident  This  function  requires  a 
minimum  of  three  valid  values  to  compute. 


SORT  Square  root  function 


CO  TY«SQRr(TX) 

Compute  the  square  root  of  TX  and  store  in  TY.  TX  and  TY  can  be  the 
same  variable.  TX  can  be  scalar,  time  series  or  paired  data,  if  a  value  of  TX 
is  undefined  resulting  values  of  TY  are  undefined.  If  a  concurrent  IF 
coTKiition  is  not  satisfied,  TY  is  unchanged.  If  a  value  in  TX  is  negative,  the 
value  in  TY  will  be  set  to  missing. 
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Name: 


Straddle  Stagger  routing  function 


Uae: 

Description: 


Name: 

Uae: 

Deecription: 


Name: 

Use: 

Description: 


SS 


CO  TY=SS(TX.NAVQ.LAG.NR) 

Route  the  uniform  time  series  variable  TX  by  the  StradcUe  stagger  hydrologic 
routing  method  and  store  It  in  variable  Vf.  Note:  Variables  TY  and  TX 
cannot  be  the  same  variable  and  variable  TX  is  permanently  modified  by 
being  lagged  by  the  SS  function.  Variable  NAVG  is  the  number  of  ordirurtes 
to  average.  Variable  LAQ  is  the  ruunber  of  ordinates  to  lag  hydrograph. 
Variable  NR  is  the  number  ol  subreaches  to  use.  The  "IP  compute  opti^ 
has  no  (  act  on  this  function. 


SSW  Willmington  District  Straddle  Stagger  routing  function 


CO  TYsSS(TXNAVG.LAG.NR) 

Route  the  uniform  time  series  variable  TX  by  the  Straddle  stagger  hydrologic 
routing  method  and  store  It  b)  variable  TY.  Note:  Variables  TY  and  TX 
carmot  be  the  same  variable.  This  function  is  similar  to  function  SS  except 
that  LAG  is  usually  zero  arxf  the  resiA  of  averaging  NAVG  values  is  stored 
in  the  NAVGed  value.  Variable  NAVG  is  the  number  of  ordinates  to  average. 
Variable  LAG  is  the  number  of  ordinates  to  lag  hydrograph.  Variable  NR  is 
the  number  of  subreaches  to  use.  For  example,  if  NAVG  is  7,  the  results  for 
ordinates  13  through  19  is  stored  hi  ordinate  19.  For  the  same  example,  if 
LAG  is  set  to  2,  the  result  is  stored  in  ordinate  21 .  The  "IP  compute  option 
has  no  effect  on  this  function. 


TAN  Tangent  trigonometric  function 


CO  TY*TAN{TX) 

Compute  the  tangent  of  TX  and  store  in  TY.  TX  and  TY  can  be  the  same 
variable.  TX  can  be  scalar,  time  series  or  paired  data.  If  a  value  of  TX  is 
undefined  resulting  values  of  TY  are  undefined.  If  a  concurrent  IF  condition 
is  not  satisfied,  TY  is  unchanged. 
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Functions 


Name: 


TS1 


Interpolatad  data  at  ragular  Intervala 


Uaa: 

Daacription: 

Name: 

Uaa: 

Description: 


Name: 

Use: 

Description: 


CO  7Y«TS1(TX.DT.TOFF) 

Derive  a  new,  regular-interval  time  series  TY  from  TX.  TY  wHi  have  a  time 
span  defined  by  the  current  time  wItkIow,  a  time  interval  of  OT,  and  a  time 
interval  offeet  of  TOFF.  TOFF  Is  time  from  the  beginning  of  the  starutard 
friterval  to  the  actual  time  of  the  data.  DT  and  TOFF  are  expressed  as  units 
of  time  as  nnT,  where  T  is  M  for  minutes,  H  for  hours,  and  D  for  days.  TX 
may  be  a  regular-  or  Irregular-Interval  time-series.  Units  and  tyiM  are 
preserved.  An  IF  oorKlition  has  no  effect 


TS2  Period  averages  at  regular  Intervals 


CO  TY*TS2(7XDT.TOFF) 

Derive  a  new,  regular-intenml  time-series  TY  from  the  instantaneous  values 
in  TX  with  a  time  span  defined  by  the  current  time  window,  a  time  intenml  of 
DT,  and  an  intenml  offset  of  TOFF.  TOFF  is  time  from  the  beginning  of  the 
interval.  OT  and  TOFF  are  expressed  as  units  of  time  as  nnT,  where  T  is  M 
for  minutes,  H  for  hours,  and  D  for  days.  TOFF  is  time  from  the  beginning 
of  the  standard  interval  to  the  actual  time  of  the  data.  TX  may  be  an  irregular 
or  regular  time-series,  but  must  be  of  type  INST-VAL  TY  will  be  typed 
PER-AVER.  Example  of  use:  deriving  daily  average  flow  from  hourly 
observed  values.  An  IF  condition  has  no  effect. 


TS3  Period  mins  or  maxs  at  regular  Intervals 


CO  TY*TS3(TX,DT,TOFF,EXT) 

Rnds  either  the  minima  or  maxima  in  TX  at  DT  intervals,  with  intenml  offsets 
at  TOFF,  and  puts  the  result  in  TY.  EXT  is  used  to  indicate  the  extreme  of 
interest:  'MIN'  or  'MAX'.  DT  arfo  TOFF  are  expressed  as  units  of  time  as 
nnT,  where  T  is  M  for  minutes,  H  for  hours,  and  D  for  days.  TOFF  is  time 
from  the  beginning  of  the  standard  irrtenral  to  the  actual  time  of  the  data. 
TY  is  typed  PER-EXTR.  Example  of  use:  finding  daily  minima  and  maxima 
from  hourly  instantaneous  observations.  An  IF  condition  has  no  effect 
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Name:  TS4  lntarp<Matad  data  at  Irregular  Intervale 


Uaa:  CO  TY=TS4(TX.T2) 

Deacriptlon:  Derh^e  a  new,  Irregular-Interval  time  series  TY  from  TX  at  the  times  for  TZOe., 

TZ  provides  the  time  pattern).  Units  and  type  are  preserved.  An  IF 
condition  has  no  effect 


Name:  TSCYCL  Tkne  aeries  cyclic  analyela 

Use:  CO  XX-TSCYCMTX,TREAT.STATFILE.DSSFILE) 

Description:  Derive  a  set  of  statistics  from  cyclic  regiriar  time  series  TX.  TX  must  be 

regular  data  at  a  IHOUR,  IDAY,  or  1MONTH  intenmi.  The  time  series  TX 
may  be,  for  example,  30  years  of  IDAY  data.  For  each  interval,  (e.g.,  1st 

day  of  the  year,  2nd  day  of  the  year . 36Sth  day  of  the  year ),  statistics 

would  be  determined.  The  14  statistics  determined  for  each  interval  are: 
maximum,  minimum,  average,  standard  deviation,  5%,  10%,  28%,  50% 
(median),  78%,  90%,  95%  percentiles,  date  of  maximum,  date  of  minimum, 
and  number  of  values  processed.  These  results  are  written  to  a  specified 
DSS  file,  and  two  text  output  flies.  The  variable  XX  is  ignored. 

TREAT  controls  the  treatment  of  newly  generated  statistic  values  when 
missing  data  occur  in  the  interval.  It  must  be  expressed  in  either  the  form 
nn#  or  nn%.  If  nn#  is  used  nn  gives  the  number  of  missing  data  values 
in  the  intenral  that  will  be  accepted  In  the  computation  of  the  new  value.  If 
nn%  is  used  nn  gives  the  percent  of  missing  data  values  in  the  interval  that 
will  be  accepted  in  the  computation  of  the  new  value,  if  the  criteria  is  met 
a  data  value  will  be  gen^ated,  if  not  a  missing  value  will  occur. 

STATFILE  is  the  base  name,  not  more  than  6  characters,  of  two  files  that  will 
be  written  containing  the  same  data  written  to  the  DSS  file.  The  first  file  wfil 
contain  ail  of  the  statistics  generated,  it  will  be  in  the  form  of  a  wide  table. 
Ks  name  will  be  'statfiie.8ts’  or  ’statfile.t’  where  ’statfiie’  is  the  portion  of  the 
name  the  user  specifies.  This  file  is  wide,  so  it  is  hard  to  view  on  some 
systems.  The  second  file  is  a  subset  of  the  first  named  '8tatfile.sum‘,  or 
'statfile.s’.  The  second  file  is  limited  to  80  columns  wide  for  easy  viewing. 
The  second  file  contains,  maximums  with  their  dates,  minimums  with  their 
dates,  average,  and  the  50%  peroentile(median)  values. 

DSSFILE  is  the  DSS  file  into  which  the  cydic  results  will  be  written.  The 
results  will  be,  for  example,  a  record  of  all  the  maximums  for  each  day  of  the 
,  year  over  the  30  year  pertod.  Fourteen  records  will  be  added  to  the  DSS  file 

specified,  each  containing  one  of  the  derived  statistics.  The  records  will  use 
the  pseudo  year  3000  for  storing  the  statistics 
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Name: 


TSHIFT 


Shift  tlma  aarfas  in  tima 


Uaa: 

Daaeription: 

Nama: 

Uaa: 

Daaeription: 


Exampla: 


COTY«TSHIFT(TX.DT) 

Dariva  a  time  aeries  TV  by  shifting  times  in  TX  by  DT.  DT  is  spedfied  in  the 
form  nnT,  where  nn  is  the  nurrfoer  of  time  units  T  and  T  may  be  M  for 
minutes,  H  for  hours,  or  D  for  days.  Units  and  type  are  preserved.  An  IF 
oondition  has  no  effect 

TTSR  Tranafbrm  tima  aariaa  to  raguiar 

CO  TY-TTSR(TXDT,T0FF,FUNCT,TREAT,TYPE) 

Derive  a  new  regular  time  series  from  existing  time-series  TX.  TX  may  be 
regular  or  irregular.  The  new  time  series  will  be  at  a  time  inten/ai  of  DT,  and 
an  intenral  offset  of  TOFF.  DT  and  TOFF  must  each  be  expressed  in  the 
form  nnT,  where  T  may  be  MIN,  HOUR,  DAY,  WEEK,  MON,  or  YEAR  ( e.g. 
IDAY).  FUNCT  is  one  of  the  following: 

INT  -  Interpolation  at  end  of  interval 

MAX  •  Maximum  over  interval 

MIN  •  Minimum  over  interval 

AVE  -  Average  over  interval 

ACC  -  Arewnuiation  over  interval 

rrO  •  Integration  over  interval 

NUM  •  Number  valid  data  over  interval 

TREAT  controls  the  treatment  of  the  new  generated  data  value  when  missing 
data  occur  in  the  interval.  It  must  be  expressed  in  either  the  form  nn#  or 
nn%.  If  nn#  is  used  nn  gives  the  number  of  missing  data  values  in  the 
intenral  that  will  be  accepted  in  the  computation  of  the  new  value.  ifnn%is 
used  nn  gives  frie  percent  of  missing  data  values  in  the  intenral  that  will  be 
accepted  in  the  computation  of  the  new  value.  If  the  criteria  is  meet  a  data 
value  will  be  generated,  if  not  a  missing  value  will  occur. 

TYPE  provides  the  user  control  to  override  how  the  interpolation  and 
processing  of  data  occurs.  TYPE  must  be  one  of  the  following:  DEFAULT, 
INST-VAL,  PER-AVEa  or  PER-CUM.  If  DEFAULT  is  used,  processing 
depends  on  the  data  type  stored  in  the  DSS  data  record.  Otherv\^ 
processing  will  be  performed  as  if  the  data  type  were  as  given  by  TYPE. 
Data  type  INST-VAL  considers  the  data  to  change  Hnearly  from  the  previous 
data  value  to  the  current  data  value  over  the  interval.  Data  type  PER-AVER 
considers  the  data  to  be  constant  at  the  current  data  value  over  the  intenral. 
Data  type  PER-CUM  considers  the  data  to  increase  from  zero  (0.0)  up  the 
current  value  over  the  interval. 

CO  TY=TTSR(TX1  DAY,0MIN>WE,1 0%,DEFAULT) 
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NaiiM: 


TTSI 


Transform  tlma  aartea  to  irragular 


Uaa: 

Daacriptlon: 

Nama: 

Uaa: 

Dascription: 


CO  TY»TTSI(TX.T2.FUNCT.TREAT.TYPE) 

Deriva  a  new  irragular  time  sarios  from  existing  time^rias  TX.  TXmayba 
regular  or  irraguiar.  TZ  is  an  existing  tima-aarias  at  the  desired  new  spacing 
ofTY.  Ait  other  parameters  are  as  shown  fc>r  the  function  TTSR. 


1/X  Invarsa  function 

CO  TY=1/X(TX) 

Divide  1  by  variable  IX  and  store  in  TY.  TX  and  TY  can  be  the  same 
variable.  TX  can  be  scalar,  time  series  or  paired  data.  If  a  value  of  TX  is 
undefined  resulting  values  of  TY  are  undefined.  If  a  concurrent  IF  condition 
is  not  satisfied,  TY  is  unchanged.  If  a  value  in  TX  is  zero,  the  value  in  TY  will 
be  set  to  missing. 
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Appendix  C  -  Programmer’s  Guide 


This  appendix  provides  histructions  on  how  a  user,  who  is  proficient  in  the  use  of  the 
Fortran  oomf^er  program  ianguage,  can  modify  existing  DSSMATH  compute  functions  and 
aeate  new  user  defined  functions.  The  DSSMATH  program  has  been  structured  to 
accommodate  the  reiatively  easy  addition  of  new  functions.  The  user  modifies  one  standard 
DSSMATH  subroutine  *USRFUN*.  in  which  a  user  defined  function  is  identified  and  then  takes 
orw  of  the  existing  function  subroutines  arxl  modifies  it  to  incorporate  the  new  function.  The 
rrtost  difficult  task  is  the  understanding  of  the  iarge  number  of  variabies  that  DSSMATH  uses  to 
store  ail  the  information  dealing  with  the  furwtion  parameters  and  data  values.  All  the  variables 
are  documented  in  the  common  blocks  in  which  they  reside. 

The  following  instructions  are  given  as  an  example  in  which  a  simple  user  defined  function 
is  created.  The  name  of  the  new  function  will  be  ABSCTX).  The  function  will  take  the  absolute 
value  of  all  values  in  the  TX  variable. 


Step  1.  •  Modify  DSSMATH  subroutine  F.USRFUN 

The  shaded  sections  of  code  reflect  the  additions  or  modifications  to  the  existing  code  to 
add  the  new  user  function  ABS.  Existing  code  is  commented  out;  if  it  needs  to  be  modified,  it 
is  duplicated  and  changes  are  made  to  the  duplicated  code. 

SUBROUTINE  USRFUN  (I8TAT) 

C  PROVIDES  A  LINK  WITH  RMCTIONS  DEVELOPED  LOCALLY 

c 

C  SUBROUTINE  AUTHOR:  DENNIS  HUFF 

C  HYDROLOGIC  ENQINEERINO  CENTER 

C  DEVELOPMENT  DATE:  8  FEB  S8 

C  REVISION  LOO: 

C  SUBROUTINE  PARAMETERS: 

C 

C  ISTAT  •  (ISTAT  -  -1)  —>  SOHETHINO  WENT  WRONG  WHILE 
C  PROCESSING  THE  COMMAND  LINE  FOR  THE  FILE  NAME. 

C  (ISTAT  -  0  )  »»>  EVERYTHING  IS  OK. 

SADD  C.IOCOM 
SADO  C.CFUNCT 
C 

C  DETERMINE  WHICH  FUNCTION  IS  BEING  USED 
C 


ELSE 


ISTAT  -  -1 

WRITE(IFOUT,*)'  ERROR  •  FUNCTION  ',CFUNEX(1:e),'  NOT  VALID’ 
ENDIF 
C 

RETURN 

END 


Appendix  C-1 
Programmers  Guide 


step  2.  -  Modify  DSSMATh  subroutine  FJVCC 

The  second  step  is  to  pick  one  of  the  existing  DSSMATH  functions  and  modify  it  to  create 
the  new  function  ABS.  in  the  exampie  beiow,  a  copy  of  the  existing  fife  F  is  modified  to 
produce  the  new  FJVBSFUN  subroutine.  F^(X}  was  chosen  because  it  is  a  function  that  is 
simiiar  to  function  ABS.  The  fbliowing  is  a  iisting  of  the  DSSMATH  function  ACC.  The  common 
biocks  have  been  turned  on  so  the  variabie  defirtitions  are  visibie. 


17 

C 

REVISION  LOO: 

ie 

C 

19 

C 

20 

C 

21 

C 

SUBROUTINE 

PARANETERS: 

22 

C 

23 

C 

. OUTPUT . 

24 

C 

25 

C 

ISTAT 

• 

(ISTAT  -  -1)  — >  SONETHINO  WENT  WRONQ  WHILE 

26 

C 

PROCEMINQ 

27 

C 

(ISTAT  •  0  )  — >  EVERYTHING  IS  OK. 

28 

C 

29 

C 

LOCAL  VARIABLES: 

30 

C 

31 

C 

NFPR 

• 

NUMBER  OF  PARAMETERS  IN  THE  FUNCTION 

32 

C 

ILQPTY 

• 

ARRAY  CONTAINING  THE  LEGITIMATE  PARAMETER  TYPES 

S3 

C 

NITSA1 

• 

NUMBER  OF  PAIRS  OF  PARAMETERS  TO  CHECK  FOR  STARTING 

34 

C 

TIMES  OF  THE  TIMES  SERIES  VARIABLES 

35 

C 

ITSA1 

• 

INTEGER  ARRAY  USED  TO  POINT  TO  WHICH  PARAMETERS 

36 

C 

ARE  TO  BE  CHECKED  AGAINST 

37 

C 

NITSA2 

• 

NUMBER  OF  PAIRS  OF  PARAMETERS  TO  CHECK  FOR  DATA 

38 

C 

TIMES  OF  THE  TIMES  SERIES  VARIABLES  OCCURING  AT 

39 

C 

THE  SAME  TIMES 

40 

C 

ITSA2 

• 

INTEGER  ARRAY  USED  TO  POINT  TO  WHICH  PARAMETERS 

41 

C 

ARE  TO  BE  CHECKED  AGAINST 

42 

c 

43 

c 
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C$AOO  C.08S0A1 
C 


THIS  COMiON  BLOCK  CONTAINS  VAAIABLES  DEALINO  WITH  THE  OSS  DATA  AND 
VARIABLE  NANES  BEING  USB)  TO  HOLD  THIS  DATA 

PARANETER  (NTS>>10.NPF-6,NSCA-18,NPFC-aO,NTSV«000,NPFV-8000, 

A  npat«nts>npf,nvm-nts>w^^nsca{ 

CONNOR  /DSSOA1/  NHEADWtNPAT) .HEA0«(10O.NPAT) ,NTYPE(NPAT) , 

A  NPNP(S.NPAT)  ,NPATH(NPAT) , ISTINE, I^TE. lETlNE, I^TE. IHORIZ(NPAT) 
A  ,TmATE<NTSV.NTS)  ,TSY(NT8V,NT8)  .PFXIN^.NPF)  ,PFY(N^,NPF) , 

A  NCURVE(NPF)  ,NiiiATA(NPAT)  ,IPATH,SO^(N^)  .IVLMiNTSV.NTS) 


INTEQER  tCAOW 

CONMON  /ROSSOA/  RBUFF 
REAL  RBUFF(NT8V*2+100) 


CONMON  /C0SS01/  VARLBL.CPNP.CPATP.CSTIHE.CSOATE.CETINE.CBMTE, 
A  CUNITS.CTYPE.PATHNH.CARYLB 
CHARACTER  VARUi(NVARl*6,CPNP(S,NPAT)*32.CPATP(S)*S2, 

A  CSTIME*4,CSOATE^7.^IMEM,C^TE*7, 

A  CUNITS(2,NPAT)*S,CTYK(2.NPAT)*S,PATHNH(NPAT)*SO, 

A  CARYLB(N^C,I^)*S 


COMMON  /LDSSOA/  LTWSET 
LOGICAL  LTWSET 

DESCRIPTION  OF  VARIABLES 


NTSV 

NPFV 

NVAR 

NHEADWO 

HEADWO 

VARLBLO 


NTYPEO 


IHORIZ 

IQFLAGO 


ISTIME 

lETIME 

ISOATE 

lEDATE 

CSTIME 

CETIME 

CSOATE 

CEDATE 

CPATPO 

IPATH 

CPNPO 


MAXIMUM  NUMBER  OF  PATHNAMES 
MAXIMUM  NUMBS)  OF  TIME  SB)IES  VECTORS 
MAXIMUM  NUMBS)  OF  PAIRED  FUNCTIONS 
MAXIMUM  NUMBS)  OF  SCALARS 

MAXIMUM  NUMBS)  OF  DEPENDENT  VECTORS  PS)  PAIRED 
FUNCTION 

MAXIMUM  SIZE  OF  TIME  SERIES  VECTOR 
MAXIMUM  SIZE  OF  PAIRED  FUNCTION  VECTOR 
MAXIMUM  NUMBS)  OF  U88)  VARIABLES 
LENGTH  OF  HEADS)  ARRAY 
OSS  DATA  HEADER  ARRAY 
USER  VARIABLE  LABELS 
FIRST  NTS  SLOTS  POINT  TO  TS  DATA 
NEXT  NPF  SLOTS  POINT  TO  PF  DATA 
LAST  NSCA  SLOTS  POINT  TO  SCALARS 
TYPE  OF  DATA  THAT  IS  HELD  BY  PATHNAME 

1  •  TIME  SERIES  REGULAR  TIME  INTERVAL 

2  •  TIME  SERIES  IRREGULAR  TIME  INTERVAL 

3  •  PAIRED  FUNCTION  DATA 

4  •  SCALAR  VARIABLE 
8  •  CHARACTER  LABEL 

VARIABLE  NUMBER  TO  APPEAR  ON  THE  HORIZONTAL  AXIS  FOR 
PLOTTING  (1  OR  2)  •  PAIRED  DATA  ONLY 
INTEGER  ARMY  TO  INDICATE  THE  FOLLOWING  FOR 
DATA  VALUES  STORED  IN  TSY()  ARRAY. 

0  •  NO  QUALITY  STATUS  IMPLIB)  > 

1  .  ESTIMATB)  VALUE 

2  -  QUESTIONED  VALUE 

3  •  MISSING  (UNDEFIIBD)  VALUE  FLAG 
TIME  WINDOW  FOR  STARTING  TIME  (MINUTES) 

TIME  WINDOW  FOR  ENDING  TIME  (MINUTES) 

TIME  WINDOW  FOR  STARTING  DATE  (DAY  COUNT) 

TIME  WINDOW  FOR  ENDING  DATE  (DAY  COUNT) 

CHARACTER  STARTING  TIME  (24  HOUR  TIME) 

CHARACTER  ENDING  TIME  (24  HOUR  TIMEJ 
CHARACTB)  STARTING  DATE  (MILITARY  STYLE) 

CHARACTER  ENDING  DATE  (MILITARY  STYLE) 


CUNITSO 

CTYPEf) 

PATHNHO 

NPATHf) 

CARYLBO 

T80ATE() 

T8Y() 


CHARACTER  ENDING  TIME  (24  HOUR  TIME) 
CHARACTB)  STARTING  DATE  (MILITARY  STYLE) 
CHARACTER  ENDING  DATE  (MILITARY  STYLE) 

LAST  PATHNAME  ENTB)ED 
POINTER  FROM  CPATP  TO  CPNP 
PATHNAME  PARTS  A-F 
LENGTH  OF  PATHNAME  PARTS  A-F 
IWITS  OF  OSS  DATA 
TYPE  OF  OSS  DATA 
FULL  PATHNAMES 
LENGTH  OF  PATHNAMES 

LABELS  FOR  PAIRED  FUNCTION  DEPENDENT  VECTORS 
TIME  SERIES  DATES  (JULIAN  DAYS  WITH  FRACTIONS) 
TIME  SB)IES  DEPEND^  VARIABLE 
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80:  C  PFX()  -  PAIRS)  FUNCTION  INDEPENOeNT  VARIABLE 

81:  C  PFY()  •  PAIRED  FUNCTION  DEPENDENT  VARIABLES 

82:  C  NCURVEO  •  NUMBER  OF  DEPENDENT  VARIABLES  PER  PAIRED  FUNCTION 

83:  C  NDATM)  -  NUHBS1  OF  ORDINATES  P61  CURVE  (PF  OR  TS) 

84:  C  LTW^  .  LOGICAL  FLAQ  TO  INDICATE  IF  THE  TIME  WINDOW  IS  SET 

SS:  C  RBUFF  •  OSS  WORK  SPACE 

S8:  C  SCALARO  •  SCALAR  VARIABLES 

S7:  C 

88:  C 

88:  C 

80:  COWnN  /TBUFFR/  XTBUFF(NTSV) ,IT8UFX(NTSV) 

01  *  c 

82:  C  XTBUFFO  •  TEMPORARY  BUFFS!  SPACE  USED  TO  STORE  A  REAL  VARIABLE 

83:  C  ITBUFXi)  •  TEMPORARY  BUFFS!  SPACE  USED  TO  STORE  AN  INTEQER  VARIABLE 

84:  C 


1 

2 


CSADO  C.IOCOM 
C 

C  . . I 

C  THIS  COMMON  BLOCK  HOLD  VARIABLES  DEALINO  WITH  INPUT  AND  OUTPUT 
C  UNIT  NUBERS  AND  OPERATING  ENVIRONMENT 

C 

COIflION  /lOCOM/  IFOUr,  INPUT,  lENVIR 
C 

C  VARIABLE  DESCRIPTION 


10: 

c 

11: 

c 

INPUT 

-  INPUT  UNIT  NUMBER 

12: 

c 

IFOUT 

•  OUTPUT  UNIT  NUMBER 

13: 

c 

lENVIR 

•  -0  •  INTERACTIVE  EXECUTION 

14: 

c 

>0  •  REAL  TIME 

15: 

c 

<0  •  BATCH 

16: 

c 

17: 

c 

1:  CSADO  C.COEPND 
2:  C 

3;  C 

4:  C  THIS  COMMON  BLOCK  CONTAINS  VARIABLES  DEALING  WITH  THE  DEPENDENT 
S:  C  VARIABU  USED  IN  THE  COMPUTE  COMMAND 
8:  C 

7:  COMMON  /COEPNO/  IC0PVA,IC0PTY,I0PVIF(5000) 

8:  COMMON  /CCOEPN/  COPVAR 

8:  CHARACTER  C0PVAR*8 

10:  C 

11:  C  DESCRIPTION  OF  VARIABLES 

12:  C  COPVAR  •  CHARACTER  VARIABLE  USED  TO  STORE  DEPENDENT  VARIABLE 

13:  C  NAME 

14:  C  ICOPVA  •  INTEGER  VARIABLE  USED  TO  SPECIFY  TWO  THINGS: 

15:  C  0  •  DEPENDENT  VARIABLE  COPVAR  HAS  NOT  BEEN  DEFINED 

16:  C  0  <  •  INDEX  OF  COPVAR  LOCATION  OF  ITS  DATA 

17:  C  ICOPTY  -  INTEGER  VARIABLE  USED  TO  SPECIFY  THE  TYPE  OF  DATA 

18:  C  STORED  IN  DEPENDENT  VARIABLE  COPVAR 

18:  C  lOPVIFO  -  INTEGER  ARRAY  USED  TO  INDICATE  IF  THE  DEPENDENT 

20:  C  VARIABLE  SHOULD  BE  PROCESSED  BASED  ON  THE  IF  FUNCTION 

21 :  C  RESULTS 

22:  C  0  •  PROCESS  THE  DATA  VALUE 

23:  C  1  .  DO  NOT  PROCESS  THE  DATA  VALUE 

24:  C  2  -  ONE  OF  THE  IF  FUNCTION  PARAMETERS  NAS  MISSING 

25:  C  AND  THEREFORE  DO  NOT  PROCESS  THE  DATA  VALUE 

26:  C 

27:  C 

28:  C 


1:  CSADO  C.CFUNCT 
2:  C 

4:  C  THIS  COMMON  BLOCK  CONTAINS  VARIABLES  DEALINO  WITH  TIB  FUNCTIONS 
5:  C  USED  IN  THE  COMPUTE  COMMAND 

6*  C 

7:  COMMON  /CFUNCT/  ICFUNP(IO) ,ICFPTY(10) ,ICFPS0(10) ,XCFUNP(10) , 

8:  .  NFUNP.LFUNOP 

8:  LOGICAL  LFUNOP 

10:  COMMON  /CCFUNC/  CFUNEX,CFUNP 

11:  '  CHARACTER*10  CFUNEX,CFUNP(10) 

12:  C  DESCRIPTION  OF  VARIABLES 

13:  C 

14:  C  CFUNEX  -  CHARACTER  VARIABLE  USED  TO  STORE  THE  FUNCTION  NAME 

15:  C  CFUW>  •  CHARACTER  VARIABLE  ARRAY  THAT  HOLDS  THE  NAMES  OF 
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16:  C 
17:  c 
18:  C 
19:  C 
20:  C 
21:  C 
22:  C 
23:  C 
24:  C 
25:  C 
26:  C 
27:  C 
28:  C 
29:  C 


FUNCTION  PARANETERS 

ICFUNP  -  1NTE0B1  ARRAY  USED  FOR  TWO  PURPOOER 

0  •  INDICATES  THAT  FUNCTION  PARAMETBt  IS  A  CONTANT 
0  <  •  INDEX  OF  FUNCTION  PARANETER  LOCATION 
ICFPTY  .  INTEOER  ARRAY  USED  TO  SPECIFY  THE  TYPE  OF  DATA  STORED 
SEE  DESCRIPTION  OF  NTYPEf)  VARIABU  IN  DSS0A1  COtflON 
ICFPSO  •  INTEOER  ARRAY  USED  TO  HOLD  SION  OF  PARAMETER  VARIABLE 
1  FOR  POSITIVE  AND  -1  FOR  NEGATIVE  PARAMETER 
LFUNOP  •  LOGICAL  VARIABLE  USED  TO  INDICATE  IF  THE  FUNCTION 

OPERATION  IS  BEING  USED  _ 

NFUNP  •  NUMBBI  OF  FUNCTION  PARAMETBIS 

XCFUNP  -  ARRAY  USED  FOR  REAL  NUMBER  CONSTANT  VALUE  OF  FUNCTION 


S6:  C 


4S:  C 

49:  C 

80: 

51: 

52: 

53:  C 

54: 

55: 

56: 

57:  C 

88: 


SET  SUBROUTINE  PERCHK  PARAMETER  CHECKS  FOR  THIS  FUNCTION 
DIMENSION  IL0PTY(1),ITSA1(1),ITSA2(1) 

DATA  NFPR.NITSAI ,N1TSA2  /I, 0,0/ 

DATA  ILOPTY  /12/ 

CALL  PERCHK(  NFPR.IL0PTY.NITSA1 .ITSA1 ,NIT8A2.ITSA2. 

.  I8TAT) 

IF(ISTAT.LT.O)  GO  TO  90C0 

CALL  8THDPV(1,ISTAT) 

GO  TO  9000 


S7:  C  PROCESS  THE  DATA 


69: 

70:  C 
71: 


DO  1000  I-1.N0ATA(ICFUNP(1)) 
CHECK  TW  IF  FUNCTION  FLAG 
IF(IOPVIF(l).IS.O) 


73: 

74:  C 

75: 

76: 

77: 

78:  C 

79: 


.  IGFLAGfl.ICOPVA}  -  0 

CHECK  TO  SEE  IF  TW  DATA  IS  MISSING 

ELSE  IF(IQFLA8(I,ICFUNP(1)).8T.2}  THEN 

.  IQFLAO(I.ICOPVA)-S 

•  T8Y(I,1c6pVA)--901 

COMPUTE  THE  ACCUMULATION 

ELSE 


.'•SsS'  s  A  *  vs*.  *•  A  *  Cs***  '  S'.v  S  *'«  jf  Y*'*'/  (  «V  S  b  •  ^  #  Cs  ^  s  . 

x'x  a  ^  )CM(  i.  -i".  o  y.'. 


82:  .  .  IQFLAa(I,IC0PVA)-0 

83:  .  ENDIF 

84:  1000  .  CONTINUE 

85:  C 

86:  N0ATA(IC0PVA)-N0ATA(ICFUNP(1)) 

07  *  0 

86:  NTYPE(ICOPVA)  -  NTYPE(ICFUNP(1)) 

89:  NHEADW(ICOPVA)  -  NHEA0W(ICFUNP(1)) 

00*  c 

91:  9000  RETURN 

92:  END 

93:  C 
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CAL-OSP.N  4/21/92 


BOOTSTRAP 
IB  BESIM 


NACM)  BEGIN 
*2 

I -FUNCTION 
I  SCREEN  START 
ENOMACRO 


MACRO  OPT  SFILE 
IPAGE 

Ua  COPY  SFILE. JOB  NEXT.BAT  >  NUL 

FINISH 

ENOMACRO 


MACRO  TMP 
IPAGE 

IR  PLT-TMP 
I  SCREEN  PLOTNOLO 
IPAGE 

ISCREEN  PLOT-OPT 
ENOMACRO 


MACRO  PCP 
IPAGE 

IR  PLT-PCP 
ISCREEN  PLOTNOLO 
IPAGE 

ISCREEN  PLOT-OPT 
ENOMACRO 


MACRO  GRO 

IIF  (  ^  .EG.  'ON  '  )  THEN 
ITEACH  a  OFF 
I  ELSE 

ITEACH  B  'ON  • 
lEWIF 

ISCREEN  PLOT-OPT 
ENOMACRO 


MACRO  PLT-TMP 
RESET .A 
TI  'V 

US  Ocvit,  California 
LE,1,Maxiaua  Taaperatuns 
LE,2,MiniauB  Taaperature 
LI.,3,Ay«rage  TaaRwratura 
LL,4,Raeord  Maxiaua  Taap 
LE,5,Racord  Miniaua  Taap 
PA  /''a/'VT»P-AIR  MAX//'‘a/*f/ 

PA  /''a/'VTEMP-AIR  MIN//*a/''f/ 

PA  /*a/*b/TEMP-AIR  AVE//*a/C0MP/ 

PA  /*a/'VTEMP-AIR  MAX-MAX//''a/08S.1917-19e6/ 
PA  /'‘a/'VTEMP-AIR  MIN-MIN//'‘a/0BS.1917-1986/ 
YM,4,116,e.5,-2,A,Rac  Max. 

SU  BLN  BLN 

IIF  (  ''a  -EG.  'ON  '  )  THEN 

OG  GRaMAJOR  STYLE-OOTTEO  CROSS-ON 

lEHOIF 

OL  CU-1  TY-2  COL-RED  SHADE-ON 
OL  CU-2  TY-2  COL-BLACK  SHADE-ON 


DL  OM  TY-2  COL-CMEN  SMM)E*OFF 
OL  OM  TY-2  COL-YELLO  SMOE-OFF 
OL  CU-S  TY-2  COL-CYAN  SMADE-OFF 
THANK  T 
FLOT 
YN  OFF 


MACHO  PLT-KP 

HESET.A 

TI  *l»  *1  “N 

US  Davis.  California 

YL.I.Pracip 

YL,2,Cui.  Fracip 

LE.I.Incraaantal  Praeipitation 

LE,2,CuHmilativa  Praeipitation 

LE.S.Horaal  Precipitation 

LE,4,Parcont  of  Noraal  Pracip 

PA  /*o/*h/PNECIP-INC//*a/'^/ 

PA  /*a/‘‘b/PRECIP-CUH//*a/''f/ 

PA  /‘'a/'VPKECIP-CUM//^a/NQHNAL/ 

PA  /''a/'VPRECIP*»OHMAL//''a/OaNP/ 

SPLIT  80 

SU  TLI  8RN  SLN 

IIF  (  '‘s  .EQ.  'ON  '  }  THEN 

06  GR-NAJOK  STYLE-DOTTED  CNOSS-ON 

lENDIF 

DL  CU-1  TYPE-1  COL-SLUE  SHADE-ON 
DL  CU-2  TYPE-2  COL-aAN  SHADE-QN 
DL  CU-3  TYPE-2  COL-RED  SHADE-OFF 
DL  CU-4  TYPE-3  COL-YELLO  SHADE-OFF 
THANK  T 
PLOT 

EHDMACRO 


MACRO  EXIT 

FINISH 

ENOMACRO 


MACRO  SETT 
ISCREEN  TIME 
IPAGE 

I SCREEN  PLOT -OPT 
ENOMACRO 


MACRO  1 

■TEACH  U  010CT1991 
■TEACH  V  T 
■SCREEN  TINE 
ENOMACRO 


MACRO  2 

■  TEACH  u  01JAN1991 
■TEACH  V  T 
■SCREEN  TIME 
ENOMACRO 


MACRO  S 
■TEACH  I  Start 
■TEACH  E  01JUN1991 
■  TEACH  ■)  2400 
■TEACH  L  U 
■TEACH  •  09 


ISCRffM  ariNE 
ISCREEN  TINE 
ENDNAOW 


NMM  E 
ITEACN  1  End 
ITEACN  I  2400 
ITEACN  E  15JAN1992 
ITEACN  L  V 
ITEACN  •  16 
ISCNEEN  6TINE 
ISCNEEN  TIME 
EHDNACNO 


NACNO  IS 
ITEACN  I  Start 
ITEACN  h  '  ' 
ITEACN  E  T-3Q0 
ITEACN  I  u 
ITEACN  •  09 
ISCNEEN  STINE 
ISCREEN  TINE 
ENDNACRO 


MACRO  RE 
ITEACN  I  End 
ITEACN  E  T-rSO 
ITEACN  I  '  ' 
ITEACN  L  V 
ITEACN  •  16 
ISCREEN  STINE 
ISCREEN  TINE 
ENDNACRO 


MACRO  XTINE 
ISCREEN  PLOT'OPT 
ENDNACRO 


MACRO  PLOTS 
RESET. A 

TI  -h  *u  *1  *¥ 

US  Davis,  California 
LE,1,NaxlauN  Taaperatura 
LE,2,Nfnfnua  Taa^ratura 
LEfS.Avarasa  Taaparatura 
LE,6,lncra«antai  Pracip 
PA  /^a/*b/TBV-AIR  NAX/Z^a/^f/ 

PA  /''a/^EMP-AIR  MIN//*a/''f/ 

PA  /^a/^b/TEMP-AIR  AVE//''a/0aNP/ 
PA  /*a/'VPRECIP-INC//''a/'‘f/ 

SU  SLN  TRI 

IIP  (  ^  .EQ.  'CN  '  }  THEN 

DC  CR-MAJOR  STYLEaOOTTED  CROSS«ON 

lENDIF 

DL  CU-I  TV>1  COL-RED  SHAOEiON 
DL  OU-2  TY-1  COL-SLACK  SNADE-QN 
OL  CU-S  TY-1  00L-6REEN  SKAOE-OFF 
DL  CU-4  TY-2  COL-BLUE  SHAOE-OFF 
TMARK  T 
PLOT 

ISC  PLOTHOLD  , 

•  PAGE 

ISC  PLOT-OPT 
ENDNACRO 


